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Abstract:  The  propagation  of  an  elastic  pulse  through  the  simple 
geometric  forms,  rods  and  plates,  is  investigated  from  a theoretical 
point  of  view  with  the  aid  of  the  small-raotion  dynamic  elastic  equations 
of  an  isotropic,  homogeneous,  dissipationloss  solid. 

The  investigation  is  restricted  to  disturbances  which  are  initially 
plane— wave  pulses  of  dilatation,  rand  formal  solutions  are  developed  by 
Fourier  transform  methods  (symmetric  and  one  sided).  The  resulting 
formal  solutions  are  developed  into  infinite  series,  the  terms  of  which 
represent  the  total  contribution  of  wave-groups  which  can  be  associated 
with  the  paths  of  minimum  transit  time  predicted  by  the  methods  of 
geometrical  optics. 

These  paths,  and  the  associated  wave  groups,  are  found  tc  be 
characterised  by  two  integers  and  n.,  which  represent  the  number 
of  times  the  thickness  of  the  plate  (or  iiamu.er  of  the  rod)  has  been 
traversed  as  a dilatation  .1  wave  and  as  a rotational  wave  respectively, 
the  variety  in  these  hs  is  found  to  r • ii  from  conversions  of 
dilatations!  wave  energy  to  rotational  w vo  merry  it  the  free  surfaces 


When  the  Poisson  ratio,  ^ , is  zero,  this  conversion  effect  does  not 
exist  for  the  disturbances  considered,  and  ail  of  the  energy  is  carried 
by  the  direct  dilatational  wave  (n^  ■ * 0) . 

The  terms  of  these  series  are  simplified  by  contour  deformation 
methods  and  are  found  to  represent  transients  with  the  minimum  transit 
timet  nredicted  by  the  methods  of  geometrical  optics.  In  the  case  of 
the  plate,  the  simplification  enables  one  to  carry  out  finite  numerical 
integrations  in  obtaining  quantitative  data  on  the  pulse  shape  and 
amplitude  of  the  disturbances  associated  with  any  specific  values  of 

n.  and  n . 
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Various  interference  effects  are  found  between  these  wave-groups 
as  t^^and/or  the  distance  of  transmission-*  o° . It  is  shown  that 
these  wave-groups  interfere  in  such  a way  that  the  total  disturbance 
becomes  asymptotically  a plane  wave  travelling  with  the  velocity 
predicted  in  the  classical  theories  of  thin  plates  and  rods. 

A comparison  is  made  of  these  theoretical  considerations  and  the 
experiments  reported  by  Hughes,  Pondrom  and  Mims,  and  their  failure  to 
identify  any  wave  groups  other  than  those  corresponding  to  the  direct 
dilatational  and  critical  .angle  paths  (all  of  which  have  - 0)  is 
explained. 


^D.  S.  Hughes,  W.  L.  Pondrom,  and  R.  L.  Mims,  Fnys.  t<ev.  75, 
15*2  - 1556,  (IV 19) • 
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I.  INTRODUCTION 


In  & recent  paper,  Hughea,  Pondrom,  and  Kims^  have  described  a quick 

^"D.  S.  Hughea,  W.  L.  Pondro*,  and  R.  L.  Mima,  Phys.  Rev.  75. 
1552-1556,  (1949). 

and  accurate  method  of  determining  the  dynamic  elastic  constants  of 
isotropic,  homogeneous  solids  frost  the  transmission  tines  of  elastic 
pulses  through  saaplee  in  the  shape  of  right  circular  cylinders. 

In  this  method  a longitudinal  elastic  pulse  is  delivered  to  one 
end  of  the  rod  by  a pieao-eleotrlc  driver  (an  X-cut  quart*  crystal),  and 
the  arrival  of  longitudinal  pulses  at  the  opposite  end  is  detected  by  a 
pieao-eleutric  detector  (another  X-out  quart -5  crystal).  In  general, 
many  pulses  arrive  at  the  detector  for  each  pulse  delivered  by  the  driver, 
and  the  transit  time  of  each  pulse  is  determined  by  an  electronic  system 
in  which  a crystal  oscillator  acts  as  the  basic  time  standard  or  clock. 

Although  only  longitudinal  excitation  and  detection  are  employed, 
some  of  the  pulses  are  transmitted  part  of  the  way  by  rotational  waves, 
and  the  velocity  of  rotational  waves  ae  well  as  the  velocity  of  dila- 
tations! waves  may  be  determined  fr-TO  these-  transit  times.  From  these 
velocities  and  the  density  of  the  material,  ail  e*“  lus  elastic  constants 
may  be  determined. 

A simple  theory  was  devised  with  the  methods  of  geometrical  optics, 
which  had  been  successfully  applied  to  many  similar  problems  in  seis- 
mology . 

Referring  to  Fig.  1,  the  pulse  delivered  by  the  driving  crystal 


1.1 


1.2 


would  be  expected  to  generate  a group  of  plane  waves  of  dilatation 

having  a continuous  distribution  of  frequencies  and  traveling  very  nearly 

parallel  to  the  free  cylindrical  wall  of  the  rod.  Part  of  the  energy  of 

these  waves  may  be  expected  to  be  converted  into  rotational  waves  at  the 

2 

free  cylindrical  boundary  since  dilatations!  waves  alone  cannot  in  general 


^H.  Poincare,  Lecons  sur  la  theorie  de  1' elasticite,  (Paris  1892), 
p.  12tff. 

Handbuch  der  Physik,  Bd.  VI  (Verlag,  Julius  Springer,  Berlin,  1928), 
pp.  323-321. 

See  also  reference  7. 


satisfy  the  free  boundary  conditions.  On  the  other  hand,  a portion  of 
the  energy  of  these  dilatationai  waves  will  reach  the  opposite  end  of 
trie  ruu  without  modification.  In  terms  of  ray  paths,  this  would  corres- 
pond to  T'.rh  rath  b~  or  a<:.  If  L is  the  length  of  the  rod  and  a 
is  the  ve  . . fy  of  dila*. ational  waves  in  the  rod,  tru  pulse  transmitted 


al onr  sutn 

t,  given  by 


arrive  at  the  detector  with  a 1 rummies  ion  tire. 


L 

1 - — 

a 


(1  = 1) 


.Returning  to  the  rotational  waves  obtained  by  transformation  at  the 
free  boundary,  such  waves  must  follow  Snell's  law  as  in  optics.  If  u is 
the  velocity  of  rotational  waves  in  the  rod,  tilts  require**  that. 

°B  . b 

sin  a,.  a • 

where  0^  is  the  angle  between  the  direction  of  travel  and  the  normal  to 

the  boundary  surface , and  0 Is  the  angle  between  hr.  r notion 

of  the  rotation  1 wave,  obtained  by  or  msfoirration  at  tdie  boundary,  ana 


the  normal  to  the  oonnorrv  surface. 
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Since  b <a,  and  in  the  present  case  ■ 90°,  the  rotational  waves 
must  travel  approximately  at  the  critical  angle  0^,  relative  to  the  normal 
of  the  boundary  surface  where  9 is  given  by 

sin  0£  - b/a  (1.3) 

Such  rotational  waves,  obtained  by  transformation  at  the  upper 
boundary,  would  eventually  strike  the  lower  boundary  and  there  give  rise 
to  a reflected  rotational  wave  and  a dilatational  wave.  In  satisfying 
Snell's  law,  these  dilatational  waves  obtained  by  transformation  at  the 
lower  boundary  must  travel  very  nearly  parallel  to  the  boundary,  and  part 
of  their  energy  would  eventually  roach  the  detector  end  of  thr  rod.  This 
energy  would  then  have  been  carried  over  such  a path  as  bqrd  where  bq  and 
rd  were  traversed  as  dilatational  waves,  and  qr  was  traversed  as  a rota- 
tional wave. 

If  D is  the  diameter  of  the  rod,  it  is  obvious  that  the  length  of 
the  path  qr  is  D esc  0,,  while  the  dis'ance  advanced  along  the  length  of 
the  rod  is  D tan  9., . Consequently  the  time  delay  in  transmission,  at, 
incurred  by  taking  the  path  bqrd  instead  of  be  is  given  by 

At  = (D/b)  esc  u.„  - (D/a)  tan  , =»  D(b-^  - (1-1) 

where  the  last  form  is  obtained  by  elimination  of  9 with  elation  (1.3). 

It  is  apparent  that  such  transformations  can,  in  general;  take  place 
a number  of  times.  Consequently,  such  paths  a3  bqrstc  and  bqrsturd  are 
possibly  end  cc.  respond  to  delays  in  transmission  of  2At  and  3At  re- 
spectively. and  each  gives  rise  to  a longitudinal  pulse  at  the  detector 
end  of  the  rod. 

The  process  of  reflection  at  the  ends  of  the  rod  can  also  rive  rise 
to  delayed  arrive  in.  in  particular,  part  of  the  energy  of  the  original 
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vlilr.tatior.fll  wave  group  can  bo  rof  nct*d  without  change  of  mole  of  trans- 
mi  salon,  once  at  each  end  of  the  rod,  and  thus  travel  the  length  of  the 
rod  three  tims3.  dimlltrly,  LL  oou.I1  Le  reflected  twice  at  each  on  i and 
travel  the  length  of  the  rod  five  times  etc. 

As  a result  of  these  two  processes,  transit  times  corresponding  to 
the  equation 

t * m(L/n)  ♦ nD(b~^  - a (1.5) 

are  to  be  expected  where  m is  any  positive  odd  integer  and  n is  any 
positive  Integer  or  zero.  Here  m is  the  number  of  times  the  length  of 
the  rod  is  traversed,  and  n is  the  number  of  delays  incurred  by  the  trans- 
fer of  mode  of  transmission  process.  Since  during  each  such  delay  the 
delayed  disturbance  moves  the  distance  D tan  along  the  length  of  the 
rod,  it  is  obvious  that  the  integers  m and  n must  satisfy  the  restri .on 

mL  nD  tan  (1  6) 

The  experimental  r<  *ults  are  in  excellent  agreement  with  editions 

(1.5)  and  (1.6),  Each  of  the  multiplicity  of  detected  longitudinal 
pulses  obtained  for  each  driving  pulse  corresponds  to  a particular  ad- 
missible combination  of  the  integers  m and  n,  and  the  appropriate  vari- 
ations are  observed  when  the  diamoter  or  length  of  the  rods  is  altered  ^ 
However,  the  relative  amplitudes  of  the  various  delayed  pulses  are  functions 
of  L,  D,  m,  n,  a,  b,  and  the  tendency  of  the  material  in  the  rod  to  dis- 
perse the  energy  in  such  a disturbance  into  random  elastic  disturbances 
(heat  motions),  and  many  of  the  arrivals  predicted  oy  equations  l 5J  and 

(1.6)  are  too  feeble  to  be  detected  under  some  experimental  conditions 
The  exact  wave  form  and  duration  of  the  driving  pulse  are  obviously 

of  considerable  importance.  The  duration  should  be  small  compared  to  each 
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of  the  times  (L/ft)  and  D(b'"'"  - a**")  in  order  thnt  the  detecta  1 pulses 

he  easily  r»*solve  I in  terms  of  arrival  times,  and  the  wave  form  should 
be  simple  but  easily  recognized  against  background  disturbances.  But, 
aside  from  these  considerations,  there  appears  to  be  no  additional  in- 
formation obtainable  from  geometric  optical  methods. 

In  an  effort  to  obtain  more  detailed  information,  the  literature 
was  searched  for  treatments  of  the  transmission  of  pulses  through  a 
finite  circular  rod  made  of  dissipationless,  homogeneous,  isotropic 
material.  This  search  was  entirely  without  success;  there  is  not  even 
an  exact  treatment  of  the  free  vibrations  of  a finite  circular  cylinder. 

■j 

. Pochharamer'5  has  given  a reasonably  complete  treatment  of  some  of  the 

3 

JL.  Pochhammer,  J.  f.  Math  (Crelle’s  Journal),  Bd.  81,  3<2h-336, 
(1376) . 

simpler  modes  of  transmission  of  simple  harmonic  waves  along  an  infinite 
circular  cylinder,  but  he  was  unable  to  adjust  his  wave  solutions  to  fit 
all  of  the  boundary  conditions  at  the  ends  of  a finite  circular  cylinder. 
However,  since  his  solutions  do  satisfy  the  boundary  conditions  on  the 
free  cylindrical  surface,  these  should  permit  a theoretical  treatment  of 
the  transfer  of  raode  process  dealt  with  above.  This  possibility  is 
exnlored  by  means  of  integral  transformations  in  sections  2 and  >4  of  this 
paper . 

Since  this  investigation  requires  the  use  of  Bessel  functions,  and 
it  was  recognized  from  the  start  that  some  form  of  appi-vxiisaticr.  might 
be  required,  it  was  found  desirable  to  investigate  pulse  transmission 
through  an  infinite  plate. 
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The  simple  hax*monic  modes  of  tranc  i Lon  in  an  elastic  plate  were 
investigated  by  Lord  Rayleigh^1  in  1389  and  by  H.  Lamb'1  in  1891  • De- 

\,ord  Rayleigh,  Proc.  Math.  Soc.  London,  20,  225-234*  (1889). 

^H.  Lamb,  Proc.  Math.  Soc.  London,  21,  p 85ff»  (1391). 

scriptions  of  the  various  modes  were  given,  and  the  existence  of  certain 
types  of  nodal  surfaces  was  discussed.  Although  Lamb^  reconsidered  the 

^H.  Lamb,  Proc.  Roy.  See.  London  (A),  93*  114-128,  (1916/17). 

problem  in  1916,  very  little  additional  progress  was  made.  However,  it 
was  noted  that  the  apparent  "wave  velocity"  belonging  to  these  modes 
could  be  infinite  and  was  often  greater  than  the  velocity  of  either  dila- 
tational  or  rotational  wav^c  in  an  unbounded  medium.  These  investigators 
werf>  aware  that  the  velocity  of  diiatational  waves  i3  the  maximum  velocity 
with  which  a transient  could  be  propagated  into  an  undisturbed  region  but 
offered  no  very  clear  picture  of  the  relationship  of  these  modes  of  trans- 
mission to  the  propagation  of  transients. 

7 

This  problem  was  taken  up  by  Prescott  in  1942  in  an  effort  to  throw 
7J.  Prescott,  Phil.  Mag.  Ser.  7,  33,  703-754,  (1942). 

g 

some  light  on  the  benavior  of  rods.  In  a review  of  this  article,  Bourgin 

^D.  G.  Bourgin,  Math.  Rev.,  4*  p 121,  (1943). 

suggested  that  these  modal  velocities  be  treated  n?  phase  velocities 
and  that  classical  group  velocity  theory  be  used  to  rationalise  tho 
act"al  rate  of  energy  flow. 
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These  modal  velocities  may  properly  be  referred  to  as  modal  phase 
velocities  but  are  certainly  not  phase  velocities  in  the  same  sense  as 
used  in  describing  propagation  in  an  unbounded  medium.  The  group  velocity 
calculated  from  these  by  Rayleigh's  classical  formula  gives,  as  usual, 
the  phase  velocity  of  the  modulation  envelope  for  a disturbance  which 
is  simple-harmonically  modulated  at  a modulation  frequency  small  compared 
to  the  frequency  of  the  carrier  mode.  These  group  velocities  are  not 

directly  related  to  the  fastest  transmission  of  an  abrupt  signal.  In 

9 

fact  Soramerfeld  and  Brillouin  clearly  demonstrated  in  1914  that  the 

^A.  Sommerfeld,  Ann.  der  Physik,  44,  177-202,  (1914)* 

L.  Brillouin,  ibid,  203,  240. 

fastest  transmission  of  an  abrupt  signal  takes  place  at  the  phase  velocity 
corresponding  to  infinite  frequency.  These  investigators  actually  dealt 
with  the  anamolous  dispersion  found  in  electromagnetic  wave  phenomena,  but 
the  results  are  capable  of  immediate  generalization  to  all  wave  trans- 
missions of  signals. 

In  1947  Cooper^  gave  & completely  satisfactory  demc nstration  that 
10J.  L.  B.  Cooper,  Phil.  Mag.  3er.  7,  3B,  1-22,  (1947). 


the  modes  of  simple  harmonic  propagation  in  a plate  were  in  quantitative 
agreement  with  the  expected  maximum  velocity  of  propagation,  the  velocity 
of  dilatational  waves  in  free  space.  His  method  was  similar  to  that  of 
Soramerfeld  and  Brillouin.  It  was  also  pointed  out  that  the  group  velocity 
suggestion  of  Bourgin  was  defective  in  several  respects.  In  elastic  pro- 


blems,  a single  velocity  of  energy  transfer  is  at  best  some  sort  of  an 
average,  as  there  must  be  at  least  two  rates  of  transfer  corresponding 
to  the  dilat&tional  and  rotational  methods  of  propagation.  In  the  inodes 
of  propagation  considered  there  are  two  plane  waves  of  dilatation  and  two 
plane  waves  of  rotation,  all  of  which  are  individually  traveling  in 
different  directions.  In  the  classical  group  velocity  theory  all  of  the 
plane  waves  are  traveling  in  the  same  direction.  * o mention  is  made  of 
the  fact  that  the  maximum  classical  group  valnjity  does  not  determine  the 
velocity  with  which  the  first  effects  of  a signal  are  transmitted. 

Because  of  the  complexity  of  his  formal  solution.  Cooper  did  not 
obtain  quantitative  information  on  anything  except  the  very  first  arrival 
time  for  a transient  disturbance.  In  Sectiort  3 there  is  considered  a 
particular  problem  in  pulse  transmission  through  a plate  which  closely 
resembles  in  nature  the  pulse  transmission  problem  in  the  rod  as  en- 
countered by  Hughes,  Pondrom,  and  Kims.  By  the  employment  of  a suitable 
series  of  mathematical  manipulations , it  is  possible  to  break  the  dis- 
turbance into  parts  segregated  according  to  the  nature  of  the  path 
(number  and  kind  of  reflections) . The  contribution  of  the  wave  groups 
making  up  each  part  is  in  turn  analyzed  in  terms  of  a group  transit  time, 
and  the  total  contribution  of  each  wave  group  is  then  expressed  in  terms 
of  a single  integration  which  can  be  carried  out  by  numerical  means  to 
obtain  detailed  quantitative  information  on  the  shape  and  amplitude  of 
all  the  possible  reflections  from  the  free  surfaces  bounding  the  plate. 
The  minimum  group  travel  tine  is  easily  shown  to  determine  the  beginning 
of  the  disturbance  carried  by  each  wave  group. 


An  extensive  effort  has  been  made  by  this  author  to  break  up  into 
similar  wave  groups  the  formal  solution  obtained  for  the  cylindrical  rod 
obtained  in  Section  2.  This  effort  was  only  partially  successful,  and  a 
brief  resume  is  given  in  Section  4. 

In  terns  of  physical  experiments,  only  the  rod  problem  can  be  set 
up  in  an  exact  m nr. or.  The  plate  problem  could  be  approximated  by  a rod 
of  rectangular  cross  section  with  one  side  of  the  rectangle  large  compared 
to  the  other  side  and  to  the  distance  from  driver  to  detector.  However, 
no  quartz  crystals  of  suitable  shape  are  currently  available  for  such  an 
experiments  As  a result,  the  concluding  section  us  devoted  to  a com- 
parison of  the  theoretical  solution  for  the  plate  problem  with  some 
typical  experimental  data  obtained  from  rods 


2.  A FORMAL  SOLUTION  FOR  THE  PROPAGATION  OF  AN  INITIALLY 
PLANE  DILATATION AL  PULSE  ALONG  A CYLINDRICAL  ROD  OF  INFINITE  LENGTH 


A direct  theoretical  approach  to  the  experiments  carried  out  by- 
Hughes,  Pondrom,  and  Mims  is  to  consider  the  transmission  of  an  initi- 
ally plane  dilatational  pulse  of  energy  along  ? rod  of  unif.  rm  circular 
cross-section.  Since  the  terminal  conditions  at  the  piezo-electric 
driver  and  piezo-electric  detector  are  both  partially  unknown  and  difficult 
(if  not  impossible)  to  handle  by  the  present  mathematical  methods,  we  shall 

| 

resort  to  a simplifying  assumption  as  to  the  behavior  of  the  driving 
crystal,  take  the  rod  to  be  infinite  in  length,  and  study  the  behavior 
of  the  average  normal  stress  on  a plane  section  of  the  rod,  normal  to  it9 
axis,  located  some  distance  along  the  rod  from  the  driver.  This  should 
enable  one  to  study  the  relative  amplitudes  of  the  various  reflections 
and  is  approximately  proportional  to  the  response  of  the  detecting  crystal. 

The  simplifying  assumptions  regarding  the  action  of  the  driving 
crystal  are  most  easily  understood  from  Fig.  2.  The  circular  rod  is 
actually  assumed  to  extend  to  infinity  in  length  along  the  Z axis  in 
either  direction  from  the  origin  01  coordinates,  and  the  driving  crystal 
is  thought  of  as  occupying  the  plane  section  of  the  rod  defined  by  Z » 0, 

R ~R0*  This  driving  piezo-electric  crystal  will  be  formally  replaced  by 
a uniform  surface  distribution  of  double-sources^  over  this  plane  section, 

^A.  E.  H.  Love,  A Treatise  on  The  Mathematical  Theory  of  Elasticity 
(Cambridge  192?)  183-189  and  304-30?. 
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2.2 


and  the  9trangth  of  the  double  sources  will  be  some  suitable  function  of 
time  after  t » 0 and  zero  before  t - 0.- 

X 2 

If  we  employ  the  physical  components  of  displacement  ^2 


12Ibid,  51-52,  89-91,  and  287-288. 


in  the  directions  of  increasing  R,  u,  and  Z respectively  and  the  corres- 
ponding components  and  of  the  body  forces  per  unit  of  mass 

which  act  upon  the  material  of  the  rod,  the  dynamic  elastic  equations  of 
a homogeneous,  isotropic  material  take  the  form 
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where  a is  the  speed  of  dilatatioral  waves,  b is  the  speed  of  rotational 
waves,  and  t is  the  time 

These  must  holi  at  ••'1.1  points  wilnin  the  tx dy  of  the  rod,  that  is 

Rfr.  R , if  the  rod  is  preswnod  ho  be  witJK.ot  f»  ult  or  v tl  **r  defects, 
o f ' 


2.3 


and  the  solutions  must  be  periodic  in  0 with  the  period  2rr  radians 
{i.e.  j^(R»0»Z.>t)  • £(R, 0 + 2rr,Z, t£]  in  order  to  be  single  valued* 

Further,  since  no  forces  Are  applied  to  the  outside  wall  of  the 
rod,  we  oust  have  the  boundary  conditions  that  the  stress  components 
3RR*  SR9*  vmn*,*k  surface  of  the  rod,  R - Kqo  These 

conditions  are  given  by 
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where  p is  the  density  of  the  rod  mst-sris!* 


Now  since  the  force  actions  produced  by  the  driving  crystal  are 
independent  of  0 and  have  r.o  component  in  the  0 direction,  it  is 
desirable  at  the  outset  to  restrict  this  study  to  those  motions  in  which 
• £ g ■ 0,  and  the  remaining  components  of  £ and  g are  independent 
of  0«  In  this  case  the  system  of  equations  (2,1)  reduces  to 
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and  the  boundary  condition e (2.2)  reduce  to 
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Further,  in  order  to  study  the  tr  .remission  of  pulses,  it  is  con- 
venient to  take  as  initial  conditions  the  simple  situation  in  which 
and  £2  and  their  first  partial  derivatives  xith  respect  to  the  time  s.:~r 
zero  throughout  the  rod  at  the  time  t • 0. 

If  the  component?;  £ and  the  body  force  per  unit  of  mass, 

are  taken  to  be  zero  for  t <■  0 and  are  assumed  to  be  known  after  t • 0 
throughout  the  rod,  a solution  of  equations  (2.3)  and  (2.4)  ia  readily 
obtained  by  the  application  of  Fourier  5-ntegral  transformations.  The 

solution,  subject  to  the  assumed  Initial  conditions,  is  unique,  and  an 

excellent  account  of  the  method  is  given  by  Titchmarsh.^ 


-JEo  C.  Titchraarsh,  Introduction  to  the  Theory  of  Fourier  Intermix 
(Oxford  1937),  particularly  Chapter  I. 
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We  vd.ll  employ  the  transformation  indicated  by  the  scheme 
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£*(R,Z,t)  ■ yL')  dZdt 
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where  £ is  the  vector  body  force  whose  components  are  some  known 
functions  of  R,Z,  and  t,  and  p is  the  transform  of  5°  The  components 

of  are  obviously  functions  of  R,  and  the  parameters  y,  and  s.  This 

integration  process  is  assumed  to  converge  for  real  values  of  y,  and 

any  complex  values  of  s such  that  the  imaginary  part  of  s is  greater 

than  some  positive  number  6.  In  symbolic  form  this  will  be  written 


Im(s)  > o>  0 (2„6) 

Similarly  for  the  displacement  vector  we  will  have  the 
transform,  U*,  given  by 
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£<K,7,t)  e^3t  **  r?^  dZdt.  ( 2.  7) 
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Now  eliminating  and  in  equation  (2.5)  with  equation  (2.3) 
we  find,  on  integrating  by  parts  tvdee  wit  "espect  tc  each  of  the 
variables  Z and  u,  the  result 
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where  it  has  been  Assumed  that  the  components  of  £ and  its  first  partial 
derivatives  with  respect,  to  t and  2 vanish  at  infinity  as  well  as  at 
t • 0.  Tills  additional  assumption  it-  in  essence  a boundary  condition 
to  which  one  is  led  by  the  physical  reasoning,  that  the  finite  energy, 
stored  in  this  system  by  the  body  forces  during  any  driving  oulse,  must 
ultimately  become  diffused  thro  igh  the  infinite  rod  30  as  to  become 
negligible  as  t -»  «>  or  Z -*  - 00  . Thus  the  displac  rnents,  strains,  and 
velocities  must  -n>0  as  t-t«°  or  r _ ©o  , 

Similarly,  regarding  the  system  of  equations  (?. u)  as  the  components 
of  a vector,  we  find  for  its  transfer?",  the  result 
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When  a solution  of  equation  (2.8)  for  u is  obtained  which  is 
regular  for  R *P.  Rq,  Im(s)>6  and  satisfies  the  boundary  conditions 
(2.9),  the  displacement  vector  £ can  be  found  by  the  reciprocal  trans- 
formation to  that  of  equation  (2*5)*  namely 

ooi  i-T  ^ ex? 

u(R,r,s)  E1'81  ~ rZ)dy'da  (2.10; 

-oO't-t T oo 

where 

t >6  >0. 

^ i — 

Although  equations  (2.8)  and  (2.9)  can  be  solved  for  u when  p is 
any  given  vet4 or  f. iction  of  R,y,  and  a,  we  need  only  consider  3uch  a 
•*  vector  p &3  would  correspond  to  the  particular  longitudinal  piezo- 

electric drive  given  the  rod  by  the  driving  crystal.  For  this  purpose 
we  must  take  ^ ■ 0 and  C,  - 0 when  Z / 0,  but  undefined  when  Z - 0. 

It  is  much  easier  for  one  to  make  definite  assumptions  about  p and  then 
interpret  these  physically,  than  to  deal  with  the  singular  body  force 
representing  a surface  distribution  of  sources.  Consequently,  we  will 
take 
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Thu?,  since  the  nr iy  non— zero  component  of  p Is  indoperdem  v f S, 
equation  (2.8;  has 
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However,  this  does  not  in  general  satisfy  the  boundary  conditions  given 
by  *y«t«  of  equations  (2.9) 

Since  the  syateoa  of  equations  (2=8)  and  (2,9)  are  linear  in  u,  one 
can  genersilise  the  above  solution  by  adding  as  & complementary  solution 
any  linear  combination  of  solutions  of  the  homogeneous  system  of  equations 
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The  solutions  of  this  equation  are  those  considered  by  Pochhararoer'*, 
and  of  these  only  the  two  solutions 
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and  the  linear  combinations  thereof  axe  regular  at  R * 0. 
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In  these  solutions  \JQ{z)  and  Jf  (z)  denote  Bessel  functions  of  the 
first  kind^  of  order  zero  and  one  respectively  of  the  complex  variable  z. 


^G.  H.  Watson,  A Treatise  on  the  Theory  of  Bessel  Functions 

(Cambridge  1944),  p«40. 
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k - |(?/br-  y“ | Im(k)  £ 0 

Two  other  solutions  can  be  found  by  replacing  the  Bessel  functions 
of  the  fir3t  kind  by  those  of  the  second  kind,  but  these  are  singular  at 
R ■ 0 and  must  be  excluded  since  the  displacements  must  be  regular  at 
R = 0 in  order  to  satisfy  equation  (2.?)  at  R - 0. 

Thus,  for  the  vector  body  force  given  by  equation  (2,11),  the  most 
general  solution  of  equation  (2.8),  regular  at  R - 0,  is  of  the  form 
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+ C^iy\j/0(hH)  ♦ C^k  Jo(xR) 


where  C,  and  C.,  are  arbitrary  numbers  independent  of  R. 

Substituting  this  result,  in  equation  (2.9)  and  solving  for  and 

C^,  one  finds  that  the  boundary  conditions  are  satisfied  if 
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By  eliminating  C.^  and  in  equation  (2.1?)  with  (2.18)  And  sub- 
stituting the  resulting  values  of  u^  and  into  the  reciprocal  Pourier 
transormation  equation  (2.10),  one  readily  obtains  expiescions  for  the 
displacements  and  In  this  connection,  it  is  convenient  to  write 
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(2.20) 


where  £zp  is  the  displacement  produced  by  the  particular  solution  of  (2.8) 

and  and  fzc  are  the  components  of  displacement  in  the  compl  omen*  &x*y 

solution  which  contains  the  influence  of  the  boundary  conditions.  It  is 

apparent  that  £zp  is  the  displacement  that  would  be  produced  if  the  rod 

were  unbounded  (R  - Oo  ) a 
o ' 

Upon  carrying  out  the  above  process  of  elimination,  it  is  clear  that 
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In  equation  (2.21)  the  integration  over  y 1*  readily  performed  in 
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terms  of  residues.  If  Z > 0,  we  may  by  Jordan's  lemma  consider  the 


^E.  T.  Whittaker  and  0.  N.  Watson,  A Course  of  Modem  Analysis 
(Cambridge  19 40),  p.  115, 


contour  of  integration  closed  by  an  infinite  semicircle  in  the  upper 
half  of  the  y-plane,  and  the  value  of  the  integral  is  the  sum  of  the 
residues  at  the  poles  in  the  upper  half  plane.  If  equation  (2.21)  is 
written  a* 
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the  evaluation  by  the  residue  at  y - — gives 
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Similarly,  we  find 
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Now,  since  the  tensile  stress  along  the  axis  of  the  rod  Szz  is 


given  by 


^ . (A  2b2)A  . 2b2  ^ 


(2.23) 


we  find,  by  assuming  the  validity  of  differentiation  under  the  integral 
sign,  that  the  ter.oile  stress  produced  by  the  particular  solution,  S _p, 

Lt  4*1 

is  given  by 
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Now,  taking 


F(t)  €i9t  dt 
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where  F(t)  - 0 for  t •<  0,  we  have  the  usual  result  for  one-sided  Fourier 
transforms^, 
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Thus,  th«  primary  disturbance  generated  by  the  singular  body  forces 
is  & plane  wave  of  dilatation  propagated  along  the  Z axis  with  speed  a 
and  a shape  determined  by  the  functional  form  of  F(t)  which  we  will  leave 
unspecified  for  the  presentc 

The  contribution  of  the  complementary  solution  to  “ay  likewise 
be  computed  by  substituting  and  into  equation  (2.23)  and  differ- 
entiating under  the  integral  signs.  This  contribution  S^.,  la  a function 
of  R,  and  the  receiving  crystal  will  respond  only  to  its  average  value « 
Since  all  displacements  are  inuspendent  of  0,  we  find,  by  averaging  S„„ 
over  a plane  section  normal  to  the  axis  of  the  rod,  the  result 
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Substituting  tho  values  of  £ R arid  from  equation  (2.22)  and 
performing  the  integration  over  ft  before  thoae  ov.  >-  y we  obtain 
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Sine*  is  independent  of  R,  it  is  obvious  that  tha  averags  value  of 

tha  antira  stress,  -S^,  gi*«n  by 


SZZ  " SZZP  * SZZC  (2o29) 

Prom  equations  (2.22)  and  (2.26),  it  is  obvious  that  tha  comple- 

2 2 

teant ary  solution  vanishes  if  a • 2b  . This  situation  corresponds  to  a 
zero  Poisson  ratio,  or  f since 

(1  - 2b2/ a2)  - — (2.30) 

1 - a 


An  actual  material  having  <r  **  0 wou.-d  be  highly  unusual  and  have 

interesting  uses.  It  would  not  violate  the  conditions  of  physical  sta- 
ir? 

bility.  Rayleigh  ' has  pointed  out  that  an  isotropic  material  is  stable 
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Lord  Rayleigh,  Proc.  Math.  Soc.  of  London,  17,  U-11,  (188?) 


if  its  Poisson  ratio  lies  ir.  the  range  1 /2.x*  <r  -1.  When  cr  - 0,  it  is 
obvious  that  no  reflections  of  the  primary  dilatational  disturbance 
considered  hers  are  produced  by  the  surrounding  walle.  Such  a material 
would  be  extremely  useful  for  making  certain  typos  of  solid  acoustic 

delay  lines o 
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The  evaluation  of  such  integrals  as  occur  in  equations  (2.22)  and 
(2. 2d)  can  be  carried  out  by  a number  of  different  procedures . Since 
the  integrands  are  even  functions  of  h and  k,  they  are  single  valued 
functions  of  s and  y.  In  addition,  they  are  regular  for  Im(»)  > 6 > 0, 

0 ^ R R_,  and  t<  0.  Thus,  the  integration  over  * is  sero  under  these 
conditions.  When  t > 0,  0^  R 4 Rq>  the  integrands  are  regular  for 
Ia(s)^  T except  for  poles  which  are  all  located  on  the  real  axis  of  the 
s piano  when  y is  real.  In  this  caBe  the  integration  over  s is  obtained 
aa  the  oum  of  the  residues  at  these  poles.  3ince  the  equation 


^(h,k,Ro)  - 0 (2.31) 

3 

ia  actually  the  frequency  wave-nussber  condition  obtained  by  Pochhamer 
for  the  normal  modes  Independent  of  0 for  the  infinite  rod,  Jt  follows 
that  this  residue  evaluation  will  express  the  integral  as  a series  of 
these  normal  modes.  Each  of  these  modes  must  then  be  integrated  over  y. 
These  processes  will  converge  very  slowly  and  are  incapable  of  showing 
the  number  of  reflections  which  have  given  rise  to  the  energy  arriving 
at  any  given  time  and  place o 

A second  approach  ia  to  expand  the  integrand  in  a series,  convergent 
for  Im(s) 6,  and  integrate  term  by  term.  This  possibility  is  a conse- 
quence of  the  relations 


Im(h)  > 7 L*(  a) 
ia(k)  > ~ Ia(s) 


(2.32) 


which  are  readily  obtained  from  equation  (2.16)  and  are  valid  for  all 


2.16 


real  values  of  y.  in  this  process,  a type  of  series  can  be  found  in  which 
the  various  type3  of  geometrical  reflected  paths  appear  as  separate  terms 
distinguished  by  appropriate  changes  of  phase  according  to  the  increased 
length  of  path*  This  method  is  pursued  with  some  success  in  Section  U,  but 
the  method  is  much  more  successful  in  lection  3 in  a similar  problem 
involving  pulse  transmission  in  a plate. 

This  problem  of  pulse  transmission  in  a flat  plate  is  much  more 
susceptible  to  formal  integration  and  is  to  some  extent  an  approximation 
to  the  present  situation.  Since  it  will  serve  as  a guide  in  the  more 
difficult  problem  of  the  rod,  it  is  considered  next. 


3o  THE  PROPAGATION  OF  AN  INITIALLY  PLANE  DILATATIONAL  PUI.JE 
THROUGH  AN  INFINITE  PLATE 


As  a parallel  investigation,  let  us  consider  an  infinite  plate  of 
uniform  thickness,  2Ro>  wnich,  if  described  in  rectangular  coordinates, 
is  bounded  by  the  two  plane  and  parallel  surfaces  X - Rq,  and  X - -R^, 
and  extends  to  infinity  in  all  directions  perpendicular  to  the  X axis.. 

We  will  assume  a uhifons  distribution  of  double  sources,  whose  strength 
depends  on  the  time  in  the  same  manner  as  in  the  previous  example,  along 
the  infinite  rectangular  strip  defined  by  Z ■ 0,  |X|  dRQ. 

The  effect  of  the  initially  plane  dilatational  pulse,  and  the  corre- 
sponding reflections  form  the  stress-free  bounding  surfaces,  car.  again  be 
readily  studied  by  calculating  the  average  stress  S,^  on  another  infinite 
strip  defined  by  Z * constant  and  |xj^  which  is  parallel  to  the  first 
strip  containing  the  double  sources  and  located  a distance,  Z,  away« 

Proceeding  in  the  same  manner  as  in  the  case  of  the  infinite  rod, 
the  dynamic  elastic  equations,  when  expressed  in  rectangular  Cartesian 
coordinates,  take  the  form 
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for  an  isotropic  and  homogeneous  solid,  whsre^,  and  are  th« 
components  of  displacement,  and  and  are  ^®  components  of  the 

body  forces  per  unit  mass  in  the  di  actions  of  increasing  X,  Y,  and  Z 
respectively. 

Since  the  planes  X • -Rq  which  bound  the  plate  must  be  stress  free, 
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SXY  * ^XZ  * ^ mus^  k®  satisfied  on  these 


surfaces.  These  mav  bo  written  as 
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where  the  symbols  p,  a,  and  b are  those  previously  defined*. 

For  the  present  circumstances,  we  wish  to  consider  the  case  in  which 

5y  ■ - 0,  and  all  of  the  other  components  of  and  C are  independent 

of  Y.  This  is  the  two-dimensional  case  discussed  at  length  by  Rayleigh^* 
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and  Lamb  ’ in  which  the  above  general  equations  reduce  to  the  system 
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with  the  boundary  conditions 
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Applying  the  Fourier  transformations  indicated  by  equations  ( 
and  (2.7)  to  this  system  with  the  assumptions  £T  ■ * 0 for  t<- 
and  fjr,  and  their  first  nartial  derivatives  are  zero  at  t ■ 0, 
obtains  the  transformed  system 
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Making  the  assumption 


PX 


0 


~2rg(3)6 

p (3.7) 
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corresponding  to  equation  (2.11)  of  the  previous  analysis,  we  find  that 
the  system  of  equations  (3 >5)  possesses  & very  simple  particular  solution 
independent  of  X which  does  not  in  general  satisfy  the  boundary  con- 
ditions (3-6). 

This  solution  is  easily  seen  to  be 

V * 0 


2yg(3) 

apsh2 


To  this  solution  we  may  add  any  linear  combination  of  solutions  of 
the  homogeneous  system  of  equations 
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(3*8) 


Since  these  equations  are  each  of  the  second  order  with  coefficients  in- 
dependent of  X,  four  linearly  independent  exponential  solutions  can  be 
found.  It  is  convenient  to  write  these  as 
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where  h and  k arc  given  by  equation  (2.16), 

Thus  we  have  as  a general  solution  of  equations  (3.5) » (3-6),  and 
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(3-10) 


where  the  four  nunberB  are  arbitrary  but  independent  of  X. 

By  substituting  this  solution  into  the  boundary  conditions  (3-6) 

and  solving  for  the  numbe-s  C , one  finds 
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Upon  substituting  the  values  of  and  thus  determined  into 
equation  (2.10),  it  is  convenient  to  write 
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where  £TD  is  again  given  by  equation  (2.21), 


t - gjl^b^a2), 
^XC  (2n)2i 


- oo  + t T 


00 


l.J  . u „ 

_ <yO  +■  i TT  - OO 


[_(k"-y2)  sln(hX) 3in(kRQ)  + 2y'c  sin(kX) sin{ hRQ)3 ay2g( s)  £ 

~(k2-y2) 2cos(hR  )sin(kR  )*4y2hk  sin(hR  )coa(kR  fl^b^sh 
_ o o o c » 


-i(st-rz) 

dyda 


and 


0 „ -2(l-2b2/a2) 

tjZC  (2n)2 


r >6  ^ o 

r0onT  f<» 
oo+t  T v! 


(3-13) 


r ; 2 "1  » -i(st-yZ) 

Uk  -r  )cos(hX)  sin(kfi  )-2h*.  c.os(kX)sLn(hR  )jay^g(s)c 


Rk2-y2)2cos(hrl  )si-'.(kR  k sin(hR  )cos(kR  j~],  b^sh 

o o o c J 


2 . 2 


(Vdn 


Here  it  is  again  found  that  the  complementary  part  of  the  solution 
and  ^ 7p  are  zero  when  the  Poisson  ratio  d-  is  zero  (a2  -•  2b2).  The 
general  similarity  to  equation  (2  22)  is  rather  obvious. 
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In  rectangular  coordinates  the  tensile  stress  component  is  given 
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Since,, in  the  present  example  the  displacements  are  indei  ?ndent  of  Y,  the 
average  value  8^  of  S ^ over  an  infinite  etrip  defined  by  Z <*  constant 
ix|;£E  R^  is  given  by 


v*- 


p(a2-  2b2)  fw  lX“Rc 

-^AiX  ♦ 


2R  az 

o 


tx  <3.: 

lo  L J 


From  this  expression  it  is  easily  shown  that 
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and  S is  given  b/  e-;untion  (7.2 7) 
ZZr 


Since  the  integrads  in  equations  (3*13)  and  (3*17)  are  even 
functions  of  h and  k,  they  are  single-valued  in  e and  y ar.d  permit  an 
•valuation  in  terms  of  the  residues  at  the  poles  which  appear  due  to 
the  vanishing  of  the  denominators.  For  t <0,  jlj  ^.Rq»  the  integrands 
are  regular  for  I (s)  >6  >0,  and  the  integral  over  a is  aero  for  all 
real  y.  Thus,  Szzc  » 0 for  t < 0.  For  t > 0,  |l|r<  R0>  the  integrands 
are  regular  for  I (a)^  V except  for  poles  which  occur  on  the  real  axle 
in  tbo  s plane  for  real  y.  The  vanishing  of  the  bracketed  part  of  the 
denominator  in  these  expressions  is  the  frequency  wave-number  condition 
for  the  simple  harmonic  modes  of  the  even  type  discovered  by  Rayleigh  \ 
and  the  evaluation  by  residues  for  t J>  0 will  result  in  an  expression 
in  terras  of  these  normal  modes.  This  process  will  be  quite  complicated 
and  is  not  capable  of  showing  what  kind  of  geometrical  oath  is  involved 
in  any  part  of  the  disturbance  arriving  at  a given  place  at  a given  time. 

This  process  pacorp.es  a practical  expedient  in  the  limiting  case 
Rq— >0,  as  all  of  the  poles  due  to  the  simple  harmonic  modes  degenerate 
into  a simple  pole.  Letting  R^->  0 in  equation  (3*17),  it  reduces  to 
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where  CQ  is  the  classical  velocit  of  longitudinal  waves  in  a thin  flat 
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E.  H.  Love,  Op.  Cit.,  p 497-498. 


This  integral  is  readily  evaluated  by  the  ssrae  steps  ao  employed  In 
connection  with  the  integral  in  equation  (2.21)  to  obtain 


Llm  S.„  - (a/C  )F(t  - \z\/C)  - F(t  - |z|/a) 
R->  0 ° 0 


which  when  combined  wit.h  equations  (2.2?)  and  (3.16)  yields 
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We  will  later  seo  that  this  is  an  asymptotic  result  valid  as  (Z/2R<) -?■  o0, 

2 2 2 

For  the  present,  let  ua  simply  note  that  for  1/2  > > -1,  a > C*  >2b 

2 2 2 

and  that  at  <7"  • 0,  a ■ CQ  » 2b  . Thus,  for  (T » 0,  the  disturbance  in 
this  limiting  case  reduces  to  the  original  dilatational  disturbance. 

For  the  practical  study  cf  the  transients  involved  in  the  comple- 
mentary solution,  a different  approach  is  much  more  successful.  By 
expressing  the  sines  and  cosines  in  equation  (3-i?)  with  exponentials 
and  multiplying  numerator  and  denominator  cf  the  integrand  by 
it  can  be  put  in  the  form 
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and  the  seriss  expansion 
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which  is  uniformly  and  absolutely  convergent  for  jwj  < 1 it  is  easily 
found  that 
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and  this  series  is  uniformly  and  absolutely  convergent  if  |w^|-*C  1 and 
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From  the  equations  (3-22)  and  the  algebraic  identity 
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which  is  valid  when  n is  a positive  integer  or  zero,  one  readily  obtains 
by  collecting  like  powers  ol  the  exponentials  the  result 
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where 
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when  n, , n„,  and  m are  positive  integers  or  iero  and  m does  not  exceed 

X A. 

the  larger  of  n^  and  n,,. 

From  the  obvious  elimination  between  equations  (3*26)  and  (3-2 4) > 
the  result 
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is  obtained.  Then  multiplying  by  the  appropriate  factors  and  collecting 
like  powers  of  the  exponentials,  we  have  the  expansion 


o o 

F 

u I 


( f J /f 


;)(n  ,-n,)l 


P 

i 

n-.' 


(r  /f  V- 
v 2-  1;  ' 


1) 


2iRo(n1h«n^k) 


(3  29) 


with  tne  stipulation  tlir. I the  leading  term  (n^-n  yO)  hus  the  value  unity. 

The  absolute  convergence  of  these  double  series  ia  easily  studied 
with  the  aid  of  a comparison  series.  It  is  obvious  that  the  terms  of 
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the  series  (3*2#)  are  in  magnitude  equal  to  or  less  than  the  corre- 
sponding terms  of  the  double  series  of  positive  terms 
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of  all  terms  in  this  series  such  that 
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nl  * n2~  n>  ^ I3  ea3Hy  shown  that 
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when  | v,  j <1  and  jv.,j  < 1.  it  is  c>ear  that  ice  series  (3.P.3)  is  abso- 
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lutely  convergent  if  j v,  j < 1 ana  |v.J<  1.  Sim*  the  smaller  magni- 

i <v 

tude  ,<f  and  v;  is  always  negative,  it  is  obvious  that  these  two  con- 
ditions are  equivalent  t • - ( l-v,  ) l ! -v  })  0 or  t.hat 
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1>|V^|[f"aR°lm(h)*  6';mc-ll7,(^  ♦ 6 "2ac^h)*Im(k^j  (3.30) 

is  a sufficient  condition  for  the  absolute  convergence  of  the  series 
( 3. 28}  and  (3-29). 

In  order  to  use  the  aeries  representation  (3-29)  over  the  range  of 
integration  involved  in  equation  (3*21)  and  integrate  term  by  term,  it  is 
desirable  that  it  be  uniformly  convergent  over  this  range  of  integration. 
It  is  easy  to  demonstrate  the  uniformity  of  its  absolute  convergence  for 
Ini(s)  ' from  the  inequality  (3-30).  It  la  evident  that  the  exponentials 
appearing  in  this  expression  have  an  appropriate  behavior  from  equation 
(2.32)  when  H 0.  but  the  quantity  (fy'f^j  requires  some  further  ex- 
amination. 

From  equation  (2.1 b)  it  is  readily  established  that 

r2  - <b\2  -a2h2)/(,2-b2>  <3-31> 


and  substituting  this  result  into  equation  (3.23)  one  finds 
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where  m^,  m and  ntj  are  the  roots  of  the  cubic  equatK 
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Investigation  of  the  roots  of  this  equation  for  finite  real  values 
of  the  Poisson  ratio  <T  reveals  that  there  is  never  more  than  one  positive 
real  root.  If  is  taken  to  be  this  root,  it  is  found  that  m^,>  1, 
EeCit^^rO,  and  tie(m^)^  0-  When  <T-  0,  all  of  the  roots  are  real,  but 

1,  one  root  is  negative,  and  the  third  is  zero. 

From  the  definitions  of  h and  k,  equation  (2.16),  it  is  readily 
established  that 
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and  since  0 and  j£  0 it  follows  that 
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Consequently,  from  the  equations  (3.32),  we  find 
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where  CR  i£>  the  speed  of  Rayleigh  surface  waves  and  is  the  positive 
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^Lord  Rayleigh,  Proc.  Lon^lc  i . Soc-,  17,  p.  4-11  (1387). 


root  of 
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Since  > 1 > b/a  for  physically  real  values  of  the  Poi3sor.  ratio,  It 
follows  that  a > b > Cu  for  physically  real  situations. 
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Collecting  the  several  inequalities,  it  is  evident  that 
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and  that  f f i thus  bounded  except  ia  y/s  -*ooor  -l/C^.  As  a 
f unction  of  y/s,  f^/f,  has  a double  pole  at  oo  and  simple  poles  at  -'1/CD. 

If  Im(s)  i 0,  it  follows  that  is  bounded  so  long  as  y is 

restricted  to  any  finite  range  of  real  values.  The  right-hand  side  of 

the  inequality  (3*30)  is  consequently  bounded  over  ell  values  of  Re(s) 

and  real  y when  Im(a)  >0  and  RQ>  0,  because  Im(h)3>  0,  IraCiO^O  and 

both  are  functions  which  vary  as  y when  y->-oo  or  when  y and  Re(s)  — > oo 

together  in  such  a way  that  y/s  ->  C and  C^l  b.  Thus,  the  right-hand 

side  of  the  inequality  (3.30)  has  a least  upper  bound  K which  is  a function 

of  Im(s)  but  uniform  over  all  values  of  Re(s)  and  real  y.  Since  Im(h) 

and  Ira(k)  are  functions  whi.ch  vary  as  Im(s),  the  right-hand  side  of  the 

inequality  ->  0 as  It.( e ) <*>  for  all  values  of  Re(s)  and  y in  the  domain 

considered,  and  tho  least  upper  bound  M ~>  0 as  lm(s)— yoo  B Consequently, 

there  exists  a positive  number  6 sucn  that  M <.  1 for  Im(s)  >6.  It  is 

thus  evident  that  the  series  (3-2 8)  and  (j.29)  are  uniformly  absolutely 

convergent  over  th*  domain  of  integration  in  equation  (3-21)  for 

Im(s)  > b > 0 The  number  b—ro*>  as  R — ^ 0 , 

o 

Substituting  the  series  (3.29)  into  equation  (3-21)  arid  integrating 
term  by  term,  one  finds 
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L->  Ls  C -R-. 

n.  O n . ~L  r; 

-L  1 ri. 


(3-37) 


3.18 


where 


o oo  t<-  T n oo 


Q (R0,z,t)  - (-i)ni 
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C ? 5 > 0 

and  the  inversion  of  the  order  of  summation  and  integration  is  justified 
by  the  uniform  convergence  of  the  integrals  involved,  as  well  as  the  series. 

Elementary  inspection  of  the  exponentials  occurring  in  these  integrals 
leads  to  the  identification  of  each  with  the  net  contribution  of  all  the 
waves  which  have  traversed  the  thickness  of  the  plate  n^+n^  times, 
times  as  a dilatational  wave  ar.d  n_>  times  as  a rotational  wave.  This 
identification  will  be  more  completely  justified  in  terms  of  the  group 
transit  time  to  be  associated  with  each  such  wave  group. 

The  integrals  in  equation  (3.38)  can  be  put  into  » much  more  useful 
.form  by  means  of  a ckrnge  of  variables  and  certain  contour  deformations 
Expressing  the  integrands  in  terms  of  the  new  variable  y instead  y by  the 
transformation  equations 

Y m r » - fjS- 

h - hs  f,  - f aU  (3.39) 
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and  the  grouo-transit  time  t.  defined  by 

n]n2 


t * yZ  + 2R  n,h  ♦ 2R  ri_k 
o J.  o 2 
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o nK.  -%ra  4 1 


C„  „ («„.z.t)  - (-D"l 


2(l-2b2/a2)2 
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O 1 _ _ r ag(s)  C _ ni  2 

(ni*n2)/  fL  * (n1-n2)/  f2j  ‘^^(^  fl^b^3a  (lrds  (3*41) 


where  1^  is  the  straight  line  (arg  y - -arg  a)  shown  in  Fig.  3-1  which 
corresponds  to  th*  real  axis  in  the  y-plane.  From  equations  (2.16), 
(3.23)  and  (3.39)  it.  is  obvious  that 

h - UA2  - r2)1/2  Re  (k)  ^ Q 

k - (1/fc2  - y2)i//2  Re  (k)  2t  ■ 0 

(3-42) 

fx  - ( k2  - y 2) 2 * 4y 2 hk 

/T“2  — 2^  2 —2  re~ 

f2  =*  (k  - y ) ♦ 4y  hk 


Next,  we  consider  trie  deformation  of  the  contour  I*'  into  a new  contour 

o 

J7  which  is  chosen  such  that  the  group-transit  time  t is  real  ard 
nin2  nln2 

l on  th**  contour.  Since  there  are  poles  of  varied 

• 1 t ranch  points  at  y *»  -(l/u)  and  combined  uoLes  and 


DO 


. vs-  a * 


3=  20 

branch  points  at  y - -(!./•),  it  will  be  convenient  to  introduce  a cut  in 
the  y plane  extending  from  y - l/a  out  the  positive  real  axis  to  oo  and 
another  cut  extending  from  y • -l/a  out  the  negative  real  axis  to  -oo, 
ill  of  the  singularities  lie  on  these  cuts  as  shown  in  Fig-  3*1  and  the 
system  of  equations  (3-4 2)  is  valid  throughout  the  plane  cut  in  this 
manner.  This  process  of  cutting  the  plane  selects  one  sheet  of  a four 
sheeted  Riemann  surface.  An  alternative  definition  of  the  sheet  selected 
is  the  conditions  Re(h)  s>  0 and  Re(k)  ^ 0. 

The  construction  of  the  required  contour  presents  no  very  serious 
problem  as  there  are  never  more  than  two  roots  of  equation  (3-40)  in  y 
consistent  with  equation  (3*42) • As  equation  (3-40)  is  an  equation  ex- 
pressing t as  a function  of  y it  is  convenient  to  define  a reciprocal 
n 
M. 

relationship  expressing  y as  a function  of  t when  t is  real . 

n « « i n H . 
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This  relationship  will  be  expressed  as  y (t)  and  will  be  defined  as 

”ln2 

the  root  lying  in  the  upper  half  of  the  y plane  when  these  two  roots  are 

complex  conjugates  and  as  the  root  with  the  smaller  magnitude  when  both 

roots  a~s  real.  This  provides  a unique  and  continuous  transformation 

from  t to  y for  real  values  of  t — 2R  Rn,/a)  + (n.,A>T]  It  is 
n.n„  n.n..  o1-  1 ?; 

J.  2 X x r _ - 

desirable  to  think  of  real  values  of  y (t)  as  lim  y (t)  + i £ 

nin2  ’ J 

because  of  the  cut-plane. 

From  the  equations  (3-42)  it  is  obvious  that  when  y is  in  the  upper 
half-plane  and  ouch  larger  in  magnitude  than  1/a  or  l/b,  h 'V  k '■'-'-iy, 
Applying  this  result  to  equation  (3-40),  it.  is  obvious  that  when  t and  y 
are  large  the  asymptotic  result 
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''J‘/[Z  - 21So(rl*"20 


(3.43) 


la  obtained,  and,  for  the  corresponding  values  of  h and  k,  the  similar 
expression 

h<vk^t/[iZ  ♦ 2Ro(ryn2)1  (3-44) 

is  obtained. 

For  large  values  of  t the  roots  of  aquation  (3.40)  are  thus  complex. 

If  we  denote  by  y (t)  the  complex  conjugate  of  y (t),  these  two 
nln2  nln2 

numbers  are  the  two  roots  of  equation  (3*40)  consistent  with  equation 

(3.42).  These  roots  decrease  in  magnitude  as  t decreases  through  real 

values  and  become  real  and  equal  at  a value  of  t which  we  will  denote  by 


t ° 
nln2 


when  the  corresponding  real  value  of  y lies  in  the  interval  -l/a 


y 2=  l/a.  Under  these  circumstances  we  may  obviously  take  1/  as  thi 

"ln2 


locus  of  y (t)  as  t decreases  from  ♦ to  t°  and  the  locue  of 


—i* 


nlr,2 


n-in? 


y (t)  as  t increases  from  t to  ♦ 00  . At  the  point  where  these 
nln2  _nln2 

tvc  loin,  it  is  obvious  that  t must  haw  a bend  point  minimun  as  y 

nln2 

varies  in  a continuous  manner.  At  this  point  we  must  have 
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Eliminating  y 


with  the  equations  (3*42),  this  takes  the  form 


2 

a 


j(n2h/  k)-*n 


(3.46) 


This  expression  is  obviously  a quartic  in  the  ratio  (h/  k)  and  can 

be  shown  to  have  only  one  real  root  (h/  k)  -i.0  and  no  pure  imaginary  roots 

for  all  real  values  of  (z/2Rq)  if  / 0,  since  a > b . Since  it  is 

easily  established  that  the  left-hand  side  of  equation  (3-46)  changes 

sign  in  the  interval  0 “^(h/  k)  -<  b^a,  it  is  apparent  that  this  positive 

root  in  (h/  k)  lie3  in  this  interval,  and  this  must  correspond  to  the 

minimum  group-transit  time  t as  it  is  the  only  root  which  satisfies 
_ nln2 

Re  (h/  k) ^ 0 which  is  valid  throughout  the  cut  plane. 

Since  it  is  easily  established  from  the  equations  (3*42)  that 

_2  b2-  &2(h/  k)2 

Y (3.47) 

aV[i-(n/k)?] 

it  is  apparent  that  for  n / 0,  the  contour  descrioed  above  crosses  the 
real  axi3  in  the  y plane  in  the  interval  -1/a  <y  -<l/a  and  thus  lies 
entirely  within  the  cut  y— plane.  Further  3tudy  reveals  that  at  this 
crossing  point  y has  the  same  sign  as  ( Z/2ft_^)  ani  that  as  (Z/2R  ) — oo , 
(h/  k)  ->0,  and  y l/a. 

The  case  n^  « 0 requires  further  study  which  is  easily  carried  out 
for  equation  (3-46)  simplifies  to 


(h/  'Y  - (b/a  Y - (1 


')  (?Y.zY 


( > . L? ) 
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Since  a > b,  this  will  have  a root  in  the  interval  b/a  > (h/Sc)  > 0 only 
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if 
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(Z/2R  ) < 5 - tan 

0 1 - (b/a)2 


(3.49) 


and  from  espiation  (3*47)  the  corresponding  value  of  y,  yQ  n is  given  by 


r°  n,  ‘ (1A)(Z/2b^0)[i  • 


(3.50) 


Substituting  this  into  equation  (3.40)  the  corresponding  minimum  group- 


transit  time  tl  is  given  by 
o J 


t° 
o n. 


UA) 


-]l/2 

2 « *2 


Z ♦ (2n  R ) 

u.  <•  0 
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(3-51) 


When  n,  « 0 and  the  inequality  (3  49)  is  eatisfied,  the  contour 
defined  above  in  terns  of  1 he  complex  roots  of  equation  (3. 40)  consistent 
with  the  equations  (3-42;  is  an  acccpta;  contour  ns  it  crosses  the  real 
axis  in  the  interval  -l/a  < Y l/a-  This  contour  is  easily  shown  to  be 
a hyperbola. 

rfhen  the  inequality  (3.49)  is  reversed,  equation  (3.48)  for  (h/k) 
hss  pure  imaginary  root3  wh-ch,  when  substituted  into  equation  (3-47); 
show  that  the  hyperbolic  contour  would  cross  the  real  axis  in  the  interval 

p 2 

l/a'^  y < 1/d  and  would  thus  cross  a cut  in  the  y plane.  This  oust  be 
avoided  by  a detour  around  the  end  of  the  cut  involved.  For  (Z/2Rq) 
positive,  this  means  a detour  around  the  point  y « l/a  and  is  readily 
accomplished . 


3.2  U 


Instead  of  crossing  the  real  axis  at  the  double  root  point  given  by 

equation  (3-50)*  t is  allowed  to  decrease  below  the  value  of  t given 

0 n2 

by  equation  (3.51)  until  the  path  cf  y (t)  now  real  hae  followed  along 

o n9 

^ ff  _ 

the  top  of  the  cut  to  the  point  l/’a,  and  the  path  of  YQ  n (0  i®  then 
followed  along  the  bottom  of  the  cut  to  return  to  the  double  root  point 
aa  t goes  through  the  same  values  in  reverse  order  The  resulting  contour 
is  sketched  iti  rig.  3>i  as  the  contour P . It  can  be  defined  as  the 

0 Ilj 

limit  of  the  contour P as  n, ->  0.  For  this  indented  contour  it  is 

nin2  1 

obvious  from  equation  (3-tO)  that  the  minimum  values  of  t and  y are  given 


by 


r,  A 2 , , 2x1/2 


t°a  n - (z/a)  v 


~'V  , ! 

y - 1/ 

o 


(3-52) 


which  is  valid  when 


(Z/2n„R  ) it  tan  9, 

O U 


The  minimum  group-transit  time  is  in  This  case  not  a bend  point  minimum. 

The  contour P may  be  thought  of  as  the  limit  of  l"1  as  n,  — > O 
oo  J r-  ir.n„  1 ^ 

L 

and  n — ^0,  and  it  obviously  runs  from  * »°  to  l/a  along  the  top  of  the 
cut  around  the  point  y *>  l/a  and  back  to  ▼ oc  along  the  bottom  of  the  cut 
when  Z ^ 0.  It  is  obvious  that  the  equations  (3-52)  hold  for  this  case 
for  all  Z ,>  0. 

The  process  of  deforming  the  original  jontour  j[^  to  the  now  contour 

P can  thus  be  accomplished  without  passing  over  any  singular  ooints  in 
Uln2 
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the  cut  y plane.  Further  It  is  possible  to  3how  from  equations  (3*36) 
and  (3*43)  that  the  integrands  vanish  as  jjy|^2nl*2n2  ^ J 

when  y — ^ ao  along  any  curve  lying  between  and  p^  n and  m is  a 
positive  number  independent  of  |y|  which  is  never  zero  except  when 

* «.  * . n tv,  no  f Vi  c ■!  nt  • crr>  ia  1 a nvar  nnv  ^ nevea  T*^  flnrt 

"1  "2  “ " 8 ±o nin2 


at  oo  are  zero,  and  the  integrals  over  y are  absolutely  convergent  at  all 
of  the  intermediate  stages  of  the  deformation  prt  ;«3a.  The  only 
difficulty  with  this  process  lies  in  letting  •.the  contours  X^D  and 
^ touch  some  of  the  singular  points.  Thcoo  difficulties  may  be 
treated  by  various  limiting  methods,  and  will  be  considered  later  as 
they  occur  in  only  a finite  number  of  terms  for  finite  values  of 
(Z/2Rq). 

Having  thus  arrived  at  the  expression 
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in  which  t-t  i3  real  at  ail  coints  in  the  range  of  integration,  the 
nln2 


order  of  integration  can  be  reversed  and  the  integration  over  s is  easily 


accomplished.  From  equation  (2.25)  and  (2.26)  it  is  readily  found  that 
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F(t)dt  - G(t) 


(3-55) 


By  applying  this  result  to  equation  (3oi+)  the  single  integral  expression 
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is  obtained.  Since  G(t)  - 0 for  t < 0,  this  expression  clearly  represents 

a transient  which  is  zero  for  t t°  , the  minimum  group-transit  time, 

nin2 

for  the  wave  groups  represented.  It  is  also  evident  tnat  the  range  of 

— —•fr 

integration  rnav  be  terminated  at  the  points  y (t)  and  y (t)  at  which 

nl  2 nln2 

t - t,  as  G(t-t  _ ) - 0 for  t <.  t. 
nl“2  ni  2 nl  2 

In  this  expression  all  of  the  waves,  which  have  traversed  the  thick- 
ness of  the  plate  n^  times  as  a diiatational  wave  and  n0  times  as  a ro- 
tational. wave,  have  been  segregated  according  to  their  transit  time.  This 
la  a highly  desirable  representation  from  a physical  point  of  view  and 
throws  some  light  on  the  more  abstruse  connections  between  wave  theory 
an  geometrical  optics. 

Considering  for  the  moment  only  real  values  of  y,  h and  It,  it  is 
evident  from  equation  (3-i*0)  and  (3-a2)  that  the  interpretation  of  these 
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variables  Ip  toms  of  geometrical  optics  is 


h - (l/a)  cos  Bp 
k - (l/b)  cos  G 

R 

y «■  (l/a)  sir.  Gn  *»  (l/b)  sin  9 


(3-57) 


where  the  notation  is  that  of  Section  1.  It  is  perfectly  obvious  that 

the  Inst  relationship  is  a restatement  of  Snell's  law,  equation  (1.2). 

A most  striking  feature  of  the  representation  in  equation  (3,56)  le 

that  the  angles  9^  and  0^  are  both  real  at  only  one  point  on  the  contour 

n } and  this  point  is  the  minimum  transit  time  point.  The  repre- 
1 2 

sentation  is  thus  rimost  enti-ely  marts  up  of  "waves"  which  have  no  geo- 
metrical meaning.  The  minimum  transit  time  wave  groups  on  the  other  hand 


have  9 and  9^  real  and  correspond  exactly  with  the  geometrical  minimum 
transit  time  paths  tc  be  expected  from  geometrical  optics 

to  demonstrate  this  equivalence,  let  us  consider  the  possible  geo- 
metric paths  in  which  the  disturbance  travels  the  thickness  of  the  plate 
times  as  a dilatations!  wave  and  times  as  a rotational  wave  obeying 
Snell's  law  where  there  is  a transfer  of  mode  and  the  law  ct  reflection 
where  there  is  no  transfer  of  mede.  Such  paths  are  shown  in  Fig.  3-2. 

It  is  obvious  that  wherever  the  disturbance  is  rotational,  it  travels 
in  a direction  making  the  angle  9^  with  the  normal  to  the  reflecting 
surfaces  of  th9  plate,  and  whorever  it  is  dilatational  with  th6  angle  0^. 

Thus,  since  the  thickness  of  the  plate  is  the  distance  moved  along 

the  plate  is  2R  tan  for  each  Uilctationai  trir>  and  2 S tan  ts  for  each 
o D o R 

rotational  trip.  The  total  distance  traveled  along  tne  plate  Z in  the 
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total  of  + n~  trip3  is  thus 


Z « 2R  n,  tan  9,.  ♦ 2R  n.  tan  9„  (3. 58) 

O J.  O O 2 K 

Making  use  of  the  transformation  equations  (3*57)  this  result  is  found  to 

be  identical  with  equation  (3*47)  which  determines  the  minimum  transit 

time  point  on  the  contour  1’ 

run,. 

± 4 

Since  the  total  distance  traveled  in  each  rotational  trip  is  2R  sec  9_. 

r o R' 

and  that  in  each  dilatations!  trip  is  2Rq  sec  0^,  the  transit  time  i3  given 
by 


2R  n.  2R  n_ 

o 1 o ? 

t * * 

n,  n~  , _ 

x 2 a cos  9^  b cos  9R 


(3-59) 


Multiplying  equation  (3- 58)  by  y and  subtracting  the  result  from  equation 
(3-59)  one  finds 


t - yZ  * 2R  n. 
n,  n , o 1 
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1 - ay  sin  9^ 

♦ 2R  n 
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1 - by  sin  0R 

1 a cos  9.. 
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which  reduces  to  equation  ( 3 * UO ) upon  using  the  transformation  equations 
\ 3 . 57 ) to  eliioinate  0^  and  0R. 

This  general  equivalence  of  the  minimum  transit  times  given  by  the 
two  methods  requires  more  study  in  the  case  “0,  for  in  this  case  the 
wave  theory  has  given  two  possible  minimum  transit  times.  The  correct, 
one  of  these  was  found  to  be  determined  according  to  whether  the  in- 
equality (3*49)  was,  or  was  not,  satisfied.  Setting  n,  * 0 in  equation 
(3. S3)  and  comparing  with  the  inequality  (3-a9),  one  finds  this  in- 
equality to  be  equivalent  to 
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(tan  eR)2^b2/(a2-  b2)  * (tan  0C)2  (3.60) 

where  0q  is  the  critical  angle  of  incidence  for  rotational  wavse  and  is 
given  by  equation  (1.3).  The  minimum  transit  time,  when  the  inequality 
(3. 49)  is  satisfied,  is  given  by  equation  (3-51)  which  is  readily  ob- 
tained setting  3 0 into  equations  (3*58)  and  (3.59)  and  eliminating 
eRi  Thus,  when  GROc,  the  minimum  transit  time  corresponds  to  rotu- 
tional  waves  travelling  a path  consisting  of  pure  reflections  a3  shown 
in  Fig.  3«3a.  Under  these  conditions,  some  interaction,  (transfer  of 
mode)  is  to  be  expected  between  the  rotational  and  dilatational  waves 
at  the  boundary  where  these  reflections  take  place.  However,  these 
correspond  to  minimum  transit  time  paths  with  n^>0,  and,  thus,  do  not 
appear  in  the  integral  expression  for  - 0. 

when  0^  > 0^,  t:  1 inequalities  (3»49)  and  (3-60)  are  not  satisfied, 
and  the  pure  reflection  path  is  not  the  quickest  path  of  transmission 
that  crosses  the  thickness  of  the  plate  n2  times  arid  aa  a rotational 
wave  each  time  In  this  case  the  critical  angle  paths,  such  as  those 
shown  in  Fig.  3 -3b,  have  a much  smaller  transit  time.  These  paths  were 
discussed  in  Section  1 in  terms  of  geometric  optics.  Setting  ra  = 1, 

L * 2,  n • n2,  and  D - 2RQinto  equation  (1.5),  one  obtains  equation  (3-52) 
which  was  derived  from  the  wave  theory  when  the  inequality  (3- 49)  was  not 
satisfied. 

It  is  obvious  that  these  two  kinds  of  path  are  identical  when 
&R  - jp.  However,  it  is  not  so  obvious  why  such  a transfer  of  mode 
process  as  trio  critical,  an ,3]  e path  is  so  important  when  the  angle  CR  for 
the  pure  ref  I eciion  pa*  h ),i -■  the  property  C»  >>  u . In  this  situation  the 
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pure  reflection  path  should  lose  no  energy  as  it  suffers  total  internal 
reflection  and  should  then  be  expected  to  make  a quite  significant  con- 
tribution to  the  overall  disturbance.  However,  in  their  experimental 
work  on  the  comparable  case  of  pulse  transmission  in  a rod,  Hughes,  Pondroa 
and  Mime  found  that  the  critical  angle  paths  explained  all  of  the  delayed 
arrivals  which  were  apparently  present. 

Geometrical  optics  offers  no  relief  from  this  apparent  contradiction, 

but  equation  (3.56)  and  the  nature  of  the  contour  P reveal  mathe- 

o 

matically  the  reason  for  thi3  experimental  result,  rfhen  these  critical 
angle  paths  give  the  minimum  transit  time,  the  contour  ^ encircles  a 
pole  at  y ■ l/a  at  the  minimum  transit  time,  and  the  sturt  of  the  corre- 
sponding disturbance  is  sharp  and  easily  identified  against  the  background. 

On  the  other  hand,  when  the  pure  reflection  path,  with  0 ©p,  gives 

K L* 

the  minimum  transit  time  or  for  the  general  case  T*  , the  contour  is 

nln2 

not  close  to  this  singularity,  and  the  disturbance  build3  up  at  a much 


alower  initial  rate  which  is  governed  by  the  values  of  dy/dt  in  the 


nln2 

neighborhood  of  the  minimum  tr.uiait.  time  point  on  the  contour,  j' 

nln2 

These  initial  rises  are  much  slower  than  those  corresponding  to  & critical 
angle  path  and  are  much  harder  to  detect  in  the  presence  of  a considerable 
background  of  other  disturbances. 

This  contrast  in  behavior  has  u distinct  .influence  on  the  use  of 
equation  (3*56)  • when  the  contours  touch  the  singui  rities  aa  in  the 
situations  where  a critical  angle  path  gives  the  first  arrival,  30019 
appear  to  analyticity  in  the  function  G(t)  is  required  so  as  to  admit 
small  contour  deformations  T.  avoid  m infinite  integrand.  This  diffi- 
culty is  easily  avoided  by  r.hc  of  a double  integral  expression. 
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Returning  to  equation  (3.5 U) > and  substituting  the  value  of  g(s) 
from  equation  (2.25) * tho  expression 
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is  obtained.  Then  using  the  familiar  result 
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the  above  expression  can  be  reduced  to  the  form 
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Since  the  functions  involved  in  this  integration  over  y are  defined  for 
all  values  of  y in  the  cut  y-  plane  except  the  actual  singular  points,  it 
is  easily  seen  that  the  contour  of  integration  may  be  any  curve  joining 

— i — # i 

the  point  y (t-t  ) to  the  point  y_  (t-t  ) which  does  not  cross  a 
“12  “l“2 

cut,  and  the  value  of  the  integral  mAy  be  found  for  all  values  of  R ,Z, 

o 

i — • . 

and  t-t  except  those  for  which  y^  n (t-t  ) falls  on  a singular  point. 
These  must  be  found  by  an  appropriate  limiting  process.  In  general,  the 
numbers  A n (R^,Z,t-t  ),  which  determine  the  response  of  the  comple- 

3-  2 , 

mentary  solution  to  a unit  impulse  delivered  at  the  time  t , can  be  found 
for  all  values  of  Rq,Z,  and  t that  are  of  physical  interest.  It  is  also 
apparent  that  they  are  capable  of  analytic  continuation  into  the  realm  of 

, i 

complex  values  of  Rq>Z,  and  ( t-t  ). 

Although  the  integral  in  equation  (3.63)  can  be  expressed  in  terras 
of  algebraic  processes  and  logarithms  for  all  of  the  integral  values  of 
n.,  and  n„,  these  expressions  are  extremely  complicated,  and  it  is  probable 
that  any  calculations  will  be  more  readily  made  by  numerical  methods  of 
integration.  As  an  example,  for  the  simplest  case  - ru  * 0,  one  finds 
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<T  is  Poissons  ratio,  auid  m^,  m ,,  and  are  the  roots  of  equation  (3-33). 
In  addition  to  the  obviously  numerous  steps  in  such  a calculation,  there 
Is  the  further  difficulty  of  having  to  compute  a small  difference  of  two 
quite  large  numbers  when  O’  is  nearly  zero.  This  iifficulty  is  present  in 
nearly  all  of  the  integral  expressions  obtained  in  this  analysis.  It  is 
due  to  the  changing  nature  of  the  singularities  located  at  y ™ -(l/a). 

As  cr->0  a pole,  located  on  another  sheet  of  the  fiiemarn  surface  as- 
sociated with  the  integrand  for  all  real  J~  / 0,  approaches  the  branch 
points  -(l/a)  as  a limit.  However,  it  is  readily  observed  that 
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avoided.  In  the  case  of  equation  (3*64)  the  obvious  approach  is  to  expand 
in  ascending  powers  of  C and  use  the  leading  term  in  this  expansion. 

The  convergence  of  the  series  (3*37)  is  obvious  for  0 and 

finite  values  of  Z and  t because  it  will  consist  of  only  a finite  mwtber 
of  non-zero  terms  under  these  conditions.  However,  the  number  of  terms 
may  be  quite  large,  and  the  individual  terms  are  difficult  to  compute.  In 
addition,  there  is  a difficulty  associated  with  the  serios  as  a whole. 

This  stems  from  the  fact  that  although  any  individual  term  has  the  property 
of  beginning  at  a certain  instant  of  time,  it  does  not  end  or,  in  general, 
even  become  small  as  t— > oo  . Indeed,  all  of  the  terms  except  that  for 
n^-n^-0  become  infinite  as  t — > co  ever,  though  the  driving  pulse,  F(t), 
has  a finite  Juration.  Thus,  the  ultimate  decay  of  the  overall  transient 
for  any  finite  value  of  Z must  be  brought  about  by  the  destructive  inter- 
ference of  the  various  wave  groups  which  have  dix'ferent  values  of  arid  rij- 
It  is  thus  necessary  to  consider  all  of  the  non-zero  terms  in  order  to  get 
an  accurate  picture  of  what  is  happening  at  any  particular  place  and  time. 

A fair  approximation  can  be  had  by  considering  only  those  wave  groups  whose 
minimum  group-transit  time  is  just  less  than  the  time  'under  consideration 
when  the  driving  pulse  is  of  snort  duration.  This  follows  from  the  be- 
havior of  the  terms  in  which  n^  ♦ is  constant.  For  these  terms  wo  find 
as  a result  of  equations  (3 <U3)  and  (3*63)  the  asymptotic  expressio 
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which  is  valid  when 


t-t'  > > (l/b)(Z2  ♦ 4n hi2)1/2 


(3.66) 


Thus,  when  t is  small  as  in  a short  driving  pulse,  we  find  that  the  sum 

of  all  of  the  terras,  for  which  n^  + n 2 is  constant,  remains  finite  as 

t — ^ OG.  Since  these  summed  groups  also  interfere  destructively  as  a 

result  of  the  factor  (-l)n  i3  equation  (3«65)»  it  thus  appears  that  a fair 

approximation  can  be  had  by  neglecting  all  terms  for  which  n^  ♦ n0  n and 

n is  such  that  the  inequality  (3-66)  is  satisfied.  However,  if  any  term 

is  neglected,  it  is  necessary  to  neglect  all  others  having  the  same  value 

of  + n2*  A somewhat  better  approximation  can  be  obtained  by  using 

equation  (3.65)  to  calculate  these  terms. 

A quite  comprehensive  interference  effect  is  also  to  be  expected  as 

Z — 7 00  , for  under  these  circumstances  the  minimum  group-transit,  time 

t approaches  the  minimum  group-transit  time  t and  the  contour 
11^2  ^ r o n0 

X*  * r as  a limit  where  n,  and  n , are  anv  finite  integers.  Thus, 

n^n?  oo  i 2 t,  r 

as  Z becomes  large,  anv  contour!1  effectively  encircles  the  singularity 

•’■n , n ,, 
i 2 
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at  y m (l/a)  for  values  of  t just  greater  than  the  minimum  value  t° 

nln2  nln2 

which  are  in  turn  just  greater  than  the  minimum  value  t^  n corresponding 

to  a critical  angle  path.  The  corresponding  contributions  to  A (R  ,Z,t) 

nln2  ° 

occur  at  times  Just  slightly  later  than  those  of  the  critical  angle  paths 

included  in  A (R  ,Z,t)  and  have  a sign  which  is  governed  by  the 
o n « o 

ni  ^ 

factor  (-1)  appearing  in  equation  (3.56).  These  successively  later 

contributions  are  qualitatively  of  about  the  same  size  as  those  of  the 

critical  angle  paths  and  as  L 00  completely  destroy  the  disturbance 

corresponding  to  the  critical  angle  path.  A similar  destruction  of  the 

primary  dilatational  disturbance  by  terms  of  the  form  h\  (R  , Z,  t)  takes 
r n^o  o 

place  as  Z — > . In  either  case,  this  destructive  interference  by  these 

slightly  delayed  but  similar  wave  trains  should  be  expected  to  ultimately 
shorten  the  wave  trains  corresponding  to  the  direct  dilatational  and 
critical  angle  paths  as  Z becomes  large.  This  is  equivalent  to  a gradual 
elimination  of  the  lower  frequency  components  in  these  wave  trains  as  Z 
increases.  In  order  to  maintain  overall  conservation  of  energy,  these 
must  reappear  in  new  wave  trains  formed  by  constructive  interference  and 
appear  at  later  times. 

The  very  complicated  undertaking  of  obtaining  an  asymptotic  series 

and  a remainder  will  not  be  attempted  here,  but  it  is  uossLbDf  to  obtain 

the  leading  term  in  such  an  expansion  in  a fairly  direct  manner.  Returning 

to  equation  (3*56)  the  integration  was  found  to  effectively  involve  only 

a finite  nart  of  the  contour  T"1  as  G(t-t  ) =*  0 for-  t t . The 

n,n_  n.n  ' n, n „ 

12  12  x 2 
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initial  and  terminal  points  were  found  to  be  y (t)  and  y (t)  res- 

nln2  nln2 

pectively,  the  values  of  y for  which  t » t. 

nln2 

For  convenience,  let  U3  consider  the  case  in  which  Z and  t oo 

in  such  a way  that  t/Z  approaches  a limit  somewhere  in  the  interval 

0 ^ (t/Z)  (l/b).  We  find  from  equation  (3«/*0)  that  y (t)  — > 

_»  nln2 

(t/Z)  4 Oi  and  y (t)  (t/2)  - Oi  for  every  finite  value  of  n,  and 

Rln2 

n^*  If  the  further  assumption  is  made  that  F(t)  - l(t)  where  l(t)  ie 
Heavisides  unit  step  function  defined  by 
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is  readily  obtained  from  equations  (3  37),  (3-54),  and  (3«5^)  in  which 
the  contour  of  integration  may  be  the  limiting  contour  or  ar.y  curve  joining 
the  initial  and  terminal  points  whior.  joen  not  cross  a cut  in  the  y -plane 
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except  at  the  singular  points.  The  double  aeries  Involved  here  ia  easily 
summed  to  yield  the  result 
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which  is  readily  integrated  to  yield 
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Since  this  is  a linear  system,  it  is  obvious  ‘ha:  he  general  result 

is 
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and  that  equation  (3«2f)  is  valid  as  Z — > ^ as  well  as  for  R — ►O:  and 

o 

thus,  it  is’  valid  as  (Z,/2R^) — ^ 00  irs  any  way  through  real  values.. 

While  the  preceding  argument  is  satisfactory  from  an  intuitive  point 
of  view,  it  is  open  to  3oine  mathematical  objection  as  a result  of  ex- 
changing the  order  of  several  limiting  processes.  In  addition,  if  the 
limiting  value  of  t/z  is  greater  than  l/b,  the  double  series  is  always 
properly  divergent  over  part  of  the  integration  contour.  A somewhat  more 
satisfactory  argument  can  be  made  directly  from  equation  (3.17)  and  will  be 


considered  in  the  next  section. 


u.  THE  TRANSFORMATION  OF  THE  FORMAL  SOLUTION  OBTAINED 


FOR  PULSE  TRANSMISSION  ALONG  A ROD 


The  successful  transformation  of  the  formal  solution  obtained  for 
pulse  transmission  along  an  infinite  plate  into  a fom  which  is  closely- 
related  to  the  considerations  of  geometrical  optics,  leads  one  to  cop**d*r 
the  possibility  of  carrying  out  a similar  transformation  upon  equation 
(2.28)  which  expresses  the  complementary  solution  for  the  similar  case  of 
pulse  transmission  along  an  infinite  rod. 

In  so  far  as  mathematical  formalities  are  concerned,  one  might  Just 
as  simply  consider  the  more  general  expression 
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where^j;(h,  k,  Rq)  is  given  by  equation  (2.19). 

This  expression  reduces  to  equation  (2.28)  when  ~j)  - l/2,  and  to 
equation  (3.17)  when  7)  «*  0.  More  generally,  when  2l)  i3  a positive 
integer,  this  expression  ia  the  solution  of  the  similar  problem  in  a hyper- 
space  having  a total  of  (2'j7  + 2)  space-like  dimensions,  and  ( 2 ♦ 1)  of 
these  snace-like  timer,  sions  making  up  the  hyp er-cv linuer  radius , 
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By  expressing  the  Bessel  functions  of  the  first  kind  in  terms  of 
those  of  the  third  kind  according  to  the  equation 

JjiU)  - a/2)  Jh”  («)  * H*  a)]  (4.2) 

where  Hj(*>  811(1  z)  are  the  two  heaeel  functions  of  the  third 

kind,  one  obtains  the  relation 
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With  the  aid  of  the  identity  (3*25)>  it  ia  easily  shown  from  equation 
(4.4)  that 
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Thus,  employing  the  aeries  expansion  (3*24),  one  obtains  the  result 
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which  is  readily  substituted  into  equation  (i«3)  and  rearranged  slightly 
to  obtain 
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with  the  stipulation  that  (h,  k,  R ) » 0 when  either  or  both  of  the 

^n,n0  o' 
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integers,  and  is  -1. 

By  using  the  3ame  methods  as  employed  in  connection  with  the  series 
(3.28),  it  is  readily  established  that  the  aeries  (4-7)  and  (4.8)  are 
absolutely  convsrgent  if 


i£yN(h,  k,  Ro)|  > \^'£\h,  k,  Ro)|  ♦ | zf*y\h,  k,  Ro)| 
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This  absolute  convergence  and  its  uniformity  are  readily  studied 

with  the  aid  of  the  asymptotic  expressions  for  Bessel  functions  of  the 
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third  kind  which  are  valid  for  large  values  of  z: 
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Since  z is  either  hP.p  or  kRy,  it  is  obvious  from  equation  (202;  that 
when  Rq)>  0,  Im(s)  is  sufficiently  large  and  positive,  and  y is  real,  the 
leading  terms  in  the  above  asymptotic  expansions  will  become  dominant  and 
the  simple  expressions 
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will  be  useful.  It  is  readily  shown  that 
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However,  an  examination  of  the  range  of  valuea  of  hRQ  and  kRQ  shows 
that  it  is  necessary  to  use  these  simple  expressions  over  the  interval 
0 < arg  z n and  a more  critical  examination  must  be  made  of  the  ex- 

VA(t)  22 

pression  for  n*)(z).  Applying  the  continuation  formulae  for  the 


22Ibid.,  p.  75- 


multiple  valued  Bessel  functions  of  the  third  kind  to  equations  (4-13) 
and  (4-14),  one  readily  finds  the  result 


2(1+M)  j cosV  tt|  € “2ljn(  z) 

3tt/2^  arg  z < tt/2 


(4.17) 


.Since  |z|  | Ira(z)  |~\  it  is  obvious  that  j£/|  and  | have  upper 

bounds  which  are  independent  of  Re(z)  and  are  monotone  decreasing  with 
the  limit  zero  as  Ira(z)  increases  through  positive  values  to  00 . Thus, 
since  a > b,  equation  (2.32)  gives  Im(z)  (R^/aJlmCs) , and  this  implies 
that  bounds  independent  of  the  real  variable  y and  Re(s)  car,  be  found 
which  are  monotone  decreasing  to  a zero  limit  as  Im(s''  increases  to  oo 
through  positive  values . Thus,  if  Im(s)  is  chosen  sufficiently  positive, 
equation  (l.ll)  is  equivalent  to  equation  (3*30),  and  the  series  (4*7) 
and  (4*3)  will  be  uniformly  and  absolutely  convergent  for  Im(s)j2L6  where 
fc  Is  positive  and <so  as  ( Vj  o°  and/or  R 


0 


It  is  unfortunate  that  the  similarity  of  the  aeries  (4*7)  and  (4.8) 
to  the  series  (3-28)  and  (3*29)  is  not  in  general  capable  of  much  greater 
extension  than  this  similarity  of  convergence.  Whereas,  each  term  of  the 
eeries  (3*28)  and  (3-29)  is  a single  valued  function  of  h and  k having  a 
finite  number  of  poles  and  an  exponential  behavior  at  infinity,  the  corre- 
sponding terms  of  the  series  (4-7)  and  (4.8)  have  a finite  number  of  poles 
and  the  same  type  of  exponential  behavior  at  infinity  only  if  -cr-Carg  (h) 
•ctt  and  -cr  < arg  (k)-Crr  and  are,  in  general,  multiple  valued  functions 
of  h and  k and  have  quite  different  characteristics  when  arg  (h)  and/or 
arg  (k)  lie  outside  the  above  open  intervals.  The  only  exceptions  to  this 
difference  in  behavior  occur  when  ^ has  such  a value  that  9in  ))n  - 0.  In 
these  exceptional  cases  the  terms  of  the  series  (4.7)  and  (4.8)  are  single 
valued  in  h and  k and  differ  in  analytical  nature  from  the  corresponding 
terms  of  the  series  (3*28)  and  (3-29)  only  in  the  number  of  poles  which 
increases  as  |])|  increases.  When  ■ 0 these  two  pairs  of  series  became 
identical. 

In  the  contour  deformations  employed  in  Section  3 in  connection  with 
the  transformation  of  equation  (3*38)  into  equation  ( 3 . 54) y arg  (h)  and 
arg  (k)  are  required  to  sweep  through  the  range  of  avlues  from 
- tan"1  [z/2(n^+n2)RQ jj  to  n ♦ tan  -1  fz/2(n^+n2)R03)  where  (z/2Rq)  is 
positive  and  the  inverse  tangent  lies  in  the  first  quadrant.  Thus,  the 
values  arg  (h)  ■ n and  arg  (k)  - tt  are  swept  over  in  the  process,  and  it 
cannot  be  employed  in  connection  with  the  term  by  term  integration  of  the 
series  (4*9)  for  sin  ^)rr  / 0 as  the  changes  in  the  behavior  for  •'-ge 
values  |hj  and  jk|  make  the  integrals  over  the  deformed  coot  .-s  oroper‘!y 
divergent. 


4 • 8 

Some  relief  from  this  difficulty  is  available  as  it  is  possible  tc 
obtain  a second  series  based  on  the  relationship 


--(1/2)  ^C2"***  H->;(z)  ♦ HyUt-2"1)] 


(4.18) 


which  ia  satisfied  by  the  Bessel  functions  of  the  first  kind  and  the 
multiple  valued  Bessel  functions  of  the  third  kind.  This  aeries  may  be 
derived  by  the  same  processes  as  the  series  (4-9)  End  will  be  written  aa 


- f f <-'>vn2wfn/  <h,  *,  .„> 

-Mh.  k,  R ) 4^4-0  nln2 


(4.19) 


n^-0  n,-0 


The  numbers  VV  (h»  k,  R ) are  easily  found  by  replacing  the 


n.n.. 
i W -‘■2 


I jw  A «■  \*/M 

functions  jjy(z)  which  occur  in  vVn  n (h,  k,  R^j , tnat  is  with  z * 


hR 


or  kRQ,  with  the  functions  £ 2nl)i  |— (z  £-2ni^  series  is  again 

uniformly  and  absolutely  convergent  for  lm(s)  > 6,  and  all  values  of  Ro(s) 
and  real  y.  However,  whereas  the  asymptotic  behavior  of  the  terms  of  the 
series  (4-9)  is  given  by 

1 


(4.20) 


* ^iO1-)  (hR^+Vn)  + n^kl^O^nyj 


4 


4-9 


for  |hRQ|  )>  ^ IVI(IVK)  -'rr-^arg  (hR^^rr 

and  jkR  I > >|V|(  1^1*1)  -n  < arg  (kR  ) < tt 


the  asymptotic  behavior  of  the  terras  of  the  aeriea  (4*19)  i»  given  by 


Vi”'™  (-1),’2Pnl4n2)/  fX  * (nl-n2)/  f2]  Pn,„2  <W 

, 0i[n1(hRo->>n)  ♦ n2(kRQ->;nQ 


(4.21) 


for  |hRol  > > |V|(|V|*1)  0 ^arg  (hRQ)  < 2n 

and  | kR  j > |V|(  |}>  (♦!)  0 < arg  (kR  )-<  2n 


Each  of  these  aeries  (4-9)  and  (4-19)  may  be  employed  over  half  of 
the  range  of  integration  to  yield 


99_  99 


SZ2C  ■ S 12  C2 

nl“0  n^-0  1 2 


(4.22) 


where  the  aeries  have  been  integrated  term  by  term  and 


1 


/*\  OH+iV 


c"* 


(R  ,Z,t)  - -i(2V*l)(l-2b2/a2)2/(2n)2 


n n ' o 
**1**2  w 


4.10 


(-D 


Vn2 


- oo 


r 


L o 


00 


2.V2.,(S)  »)/  1(.wz) 

bV*  W»l"2(  ' ' °)£  * 


(4.23) 


r°  2.Va8(.)  (h_  k_  .1(awz, 


bW 


nln2 


oO 


ds 


Z>  $ 

Upon  changing  to  tha  new  variables  y,  h,  and  k definod  by  the  system 
of  equations  (3«39)»  it  is  readily  found  th  1 tne  contours  for  the  y inte- 
gration can  be  deformed  to  yield 

\OOii~V 


n n (R.Z.t)  - -l(2>?»l)(l~2bfc/a')'/(2n)2  (-1) 


r 


n!n2  ° 


uJ  + * TT 


Vn2 


2a3  y2ag(3)\A/^)  _ -_i8(t-^Z) 

-T= 3 Wn  n <9h’sk'V  6 <V 

b4  lv  R nl  2 0 


L Jn: 


(4-24) 


2s3  y2ag(s)  . A/Wf  / _ _ 

- W „ (sh,sk,R  )6  Mt_rZ/  dy 


b^  h3R  nln2 

o 


da 


"<"R 


where  the  contour  X7  starts  at  y « 0 and  runs  along  the  positive  real 

nln2 

axis  until  it  meets  the  contour  of  fiu.  (3.1).  It-  then  goes  to  eO 

rln2  ' _ 


4.11 


along  the  branch  of  the  contour  I1  which  lies  in  the  upper  half  of 

n-jn 

— 

the  y-plane.  The  contour  JT  ia  the  mirror  image  of  r in  the  real 

nln2  nln2 

axis,  but  it  begins  at  oc  on  the  lower*  half-plane  branch  of  * and 

nln2 

snds  at  y * 0. 

Difficulties  are  again  encountered  when  n^  ■ 0 ae  the  contours  may- 
touch  the  singular  points  y • i(l/*),<  -(1/b),  and  -(l/CD)  at  which  one  or 
more  of  the  poles  of  the  integrands  in  equation  (4.24) , considered  as  a 
function  of  s,  may  recode  to  Stated  in  other  words,  when  Im(s)  is 

sufficiently  large  and  positive,  the  integrands  considered  as  • function 
of  y have  a finite  number  of  poles  (the  number  dependent  on  l)  ) clustered 
about  the  branch  points  y • -(l/a)  and  -(l/b)  and  one  in  the  Immediate 
vicinity  of  each  of  the  points  y • -(l/C^).  These  poles  approach  the 
points  mentioned  as  limits  as  Im(s)  — > 00  , and  by  choosing  6 sufficiently 
large,  the  contour  deformation  con  be  accomplished  without,  passing  over 
these  singularities  so  long  aa  / 0.  The  integrals  with  * 0 cannot 
be  handled  in  this  way  with  a finite  value  of  6.  Various  methods  of  treat- 
ing these  integrals  suggest  themselves-  The  most  obvious  of  these  is  to 
represent  these  integrals  as  the  limit  of  the  more  general  integrals  as 

nj  ->  0. 

Substituting  . e value  of  g(s)  from  equation  (2„25),  and  inverting 

the  order  of  integration } one  obtains  the  result. 


r&)  . 


(K,/VJ  « 


j 


a«  , / 's  • 

Mr  n ) F(t  )dt 


(4.25) 


4.12 


where 


Atfn  (R  - -i(2^*l)(l-2b2/a2)2/(2rr): 

• m ‘*r«  O 


f 


T*2 


<30  4-fc  T 


, -,n1*n2 

K-J-) 


W-  <jo  + i T 


L U 
r^e 


^ h3ft 


2s3  y2a 
bU  h3R 


n n 

nl"2 


d/- 


(4.26) 


VV^^sh,  3k,  R )0-i*(t-t  -rz) 


nln2 


ds 


I _ 

vanishing  of  the  integral,  over  a for  t-t  <<  t 


nln2 


and  the  finite  limits  on  the  integration  over  y are  justified  by  the 

This  vanishing  of 

the  integral  over  s in  turn  follows  from  the  fact  that  the  integrand  is 
regular  as  a function  of  a for  any  fixed  y and  Im(a)  ^6  when  t-t  t 

nln2 

The  integration  contours  for  the  variable  y can  now  be  deformed  eo  as  to 

avoid  the  necessity,  when  - 0,  of  touching  the  points  y ■ -(l/a),  l(l/b), 

and  -(l/CD)  except  when  one  of  theBe  points  must  be  the  terminal  point  on 
K 

the  contour.  In  such  a case  a limiting  process  may  still  be  required. 

Although  the  numbers  n (H0>Z,  t-t  ) are  easily  shown  to  be  zero 
for  t-t  (n^/a  + n^/b)?^,  they  do  not  appear  to  be  zero  for  t-t  ■<  n 

except  for  the  situations  when  sinj^rr  ■ 0.  In  these  special  cases  the 
two  series  (4-9)  and  (4*19)  have  identical  terms,  and  the  Bessel  functions 
can  be  expressed  in  terms  of  algebraic  functions  and  exponentials.  Thus, 
equation  (£,.26)  takes  tho  form 


A^n  OL.Z.t-t')-  -i(2V^l)(l-2b2/a2)2/(2rr)2 
nln2  ° 

T*  (ah,  sk»  R )£ 
h^R  nln2  ° 

•y.*  /t-t') 

r"iV6  J sin^n  - 0,-[>  6 


r 


00  f i "C” 


r v* 


...  T -2  -ia(  t-t-C  ) 

25'  y a -p^;  . ni«2  - 

-5 — 7T—  (ah,  sk,  R )C  4 dfde 

bU  h3R  nln2  0 


4-13 


(4«27) 


tW) 

where  / ia  a rational  function  of  ah,  ak,  and  R , and  the  contour 

nln2  _ 0 

of  integration  over  y may  be  or  any  other  curve  joining  the  point 

nln2 

y (t-t  ) to  the  point  y (t-t  ) which  does  not  cross  the  cuts  in  the 
nln2  nlr-2 

cut  y-plane  as  3hown  in  Fig.  3*1*  Since  the  integral  over  s vanishes  for 

t A(V^  , » 1 . o 

t-t  ^ t_  _ , it  is  obvi.ous  that  (R  ,Z,t-t  ) ■ 0 for  t-t  -<C  t 


*-  r> 

ln2 


n1n2  0 


nln2 


and  may  thus  be  associated  with  the  same  type  of  geometrical  path  as  the 

impulsive  response  An  n (R0,Z,t-t'  ) was  in  Section  3 whenever  3ini)rr  - 0. 

The  failure  of  the  correspondence  between  the  geometrical  paths  and 

the  impulse  responses,  A (R  ,Z,t-t  ),  when  3in  rr  jt  0,  is  rather 

n,  n0  0 

-L  *- 

disconcerting  since  the  asymptotic  properties  displayed  by  equations  (4.20) 

and  (4*21)  would  lead  one  to  oxpect  these  impulse  responses  to  be  associated 

with  waves  which  have  traveled  the  distance  2n,R  as  a dilatational  wave 

1 o 

and  2n_R  as  a rotational  wave  regardless  of  tho  value  of  1/  - 
<.  o 

The  failure  is  readily  traced  to  the  integration  over  y Li  the  interval. 

0 — y c?  y°  which  occurs  in  the  two  integrals  appearing  in  each  of  the 
nln2. 

equations  (+.24)  and  (4«2 6)  when  the  contours  J"1  and  J1"  are  employed, 

n,n 

L s*  l *4 


4*14 

When  sin>JTT  - 0,  the  contributions  from  this  part  of  the  range  of  inte- 
gration are  equal  from  the  two  integrals  but  opposite  in  sign  and  cancel 


out.  When  sin))  rr  / u,  the  contributions  are  unequal  and  do  not  cancel 


each  other. 

A study  of  the  steps  taken  in  obtaining  equations  (4*2 4)  and  (4.26) 

reveals  that  the  process  can  be  generalized  to  the  extent  that  the  contours 

F+  and  P do  not  have  to  be  terminated  at  y ■ 0,  and  could  have  been 
nln2  nln2 

terminated  at  any  other  common  point  yc  lying  in  the  open  interval  of  real 
values  -(l/a)<  yc  ■<  (l/a) . This  follows  from  the  uniform  and  absolute 
convergence  of  the  series  (4*9)  and  (4.19)  t.o  a common  Bum  when 
Im(s)>  6 and  y is  on  the  real  axis  in  the  interval  -(l/a)<  y <(l/a) . 

Terminating  the  contours  at  such  a common  point  y will  in  general  alter 

Cfc))  a A))  ? 

(R  ,Z,t)  and  A ' (R  ,Z,t-t  ) but  will  not  alter 
n,n  o ' 'n,n„  o 

(0)  1 ^ 1 ^ 

as  the  3um  of  the  changes  due  to  altering  the  contours  n will 

exactly  cancel  the  aum  of  the  changes  due  to  altering  the  contours  n 

/ 12 

for  any  value  of  l).  Thus,  in  any  calculation  of  S‘  ~ we  may  pick  y at 

any  convenient  point  on  the  above  interval.  Since  the  minimum  group-transit 

time  on  the  thus  generalized  contours  P and  T1  ” will  in  general 

nln2  nln2 

be  a function  of  y , this  choice  may  be  used  to  advantage. 

For  example  noting  that 


2Z  Dn'  (Rrt,Z.t) 


ZZC  * 

n2.“U  2 


(4^28) 
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where 


(4-29) 


and 


n^)  , 

D (R  ,z,t-t  ) 

n o 


eo 

f A« 

JLj  An  n (R0»Z>t-t  ) 
n, -0  nln2  ° 


(4-30) 


— frO  1 

we  find  on  letting  y (1/a)  - 0 that  [3  (R  »Z,t-t  ) is  0 when 
, C n2 
t-t  -4.  t^  ^ where  t^  n is  given  by  equation  (3-52).  It  is  thus  apparent 

that  although  the  correspondence  with  the  general  geometric  paths  is  in- 
complete for  sinv>rr  / 0,  it  is  alleys  possible  to  split  S^.,  up  into  parts 
which  correspond  to  direct  dilata^lonal  and  critical  angle  paths  such  as 
are  observed  for  rods. 

The  differences  in  the  behavior  of  the  series  noted  above  lead  one 
to  consider  the  asymptotic  properties  of  S„„„  as  (Z/R  ) 00  . This 

may  be  accomplished  di. ectly  from  equation  (4.1) 

4.3  in  Section  3»  it  is  convenient  to  take  F(t)  « l(t)  where  the  unit 
step-function  is  defined  by  equation  (?.6?).  The  corresponding  transform 
g(s)  is  given  by 


_\ 

6\°/ 


_ ! 


(i/3) 


(4.31) 


Substituting  this  value  into  equation  (.4.1)  and  changing  from  the 

1 1 

variables  y and  s to  tho  variables  y and  s defined  by 


y • z r 


8 - 2 8 


(4.32) 


one  finds 

oo  + i'V 

0,00  s« 

(l-2b2/»~r  I 


gC  + CT  V - & 


Ro/Z)e~^9  (t/Z)  " Y 


>if(h')3(Ro/Z)^);(h,Jk',ao/Z) 


wher» 


h - (a  /a)2  - ( y )2 
k (s  /b)2  - (r  )^j 


ni/2 


It  is  then  readily  shown  that 


L1“  3zzc  ■ 7-n 
(?./no)+»  ( ' 


o o+iT  P <7° 


I 

ir  J. 


— 11 £-iO>  (t/z)  - r J < • 

(.’)2-(^’)2  <.')2-  U„/)2  ^ 


ar» 


(4-35) 
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where  cv  is  the  positive  root  or' 


or 


a2(V  *1)  - b2(2V*l) 
a2(V*  1/2)  - 2b2 


, 2C-P  ♦ 1 

2b 

2 o~i>  ♦ (l  - «r) 


(4.36) 


(when  expressed  in  terms  of  b,  l)  , and  the  Poisson  ratio,  £T, ) and  it  is 
assumed  that  t increases  with  Z in  such  a way  that  (t/Z)  approaches  a 
finite  limit. 

The  integration  is  readily  carried  out  in  terms  of  residues  with  the 
aid  of  Jordan's  lemma  to  obtain 

Lim  sgj,  - -l(t  - |z|/a)  * (a/cy  ) l(t  - |z|/cj>) 

(Z/Ro)V«o 


Since  we  are  dealing  with  a linear  system  ar.d  such  step-functions  as 
l(t)  may  he  combined  in  a linear  fashion  to  obtain  a suitably  arbitrary 
function  F(t),  it  is  obvious  that  the  general  result  is 

Lim  g£c  - -F(t  - j Z | / a ) ♦ (a/cp)  F(t  - \z\/c^) 

(Z/Ro)->«o 

and  that 

Lim  3ZZ  - Lim  (3zzp  * S^)  - (a/c-p  ) P(t  - Ui/cy  ) (4  37 

(Z /HQy*oo  (z/Ro)-»«> 

for  any  real  value  of  C-y. 


4.18 


This  result  implies  that  at  large  values  of  (z/2Rq)  the  disturbance 

appears  to  have  travelled  the  entire  distance  Z with  the  speed  Cy . It 

is  easily  shown  from  equations  (4*36)  and  (2-30)  that  2.  2b  ^ 

for  all  real  values  of  ->)  -(1/2)  and  values  of  Poissons  ratio  (Tin  the 

range  of  physical  stability  1/2  2:  OT  £.  -1.  In  addition  for  CT  - 0 one 
° 2 2 

finds  a ® c « 2b  , and  the  complementary  solution  vanishes  identically 
as  would  be  expected. 

For  -y)  • 0,  equations  (4*36)  and  (4*37)  are  identical  with  equations 
(3-19)  and  (3*20)  respectively.  Por  V ■ 1/2  we  find 

2 2 

2 3&  - 4b  ^ 2 

(°  Vfc.)  2 -y-  - 2b4£(l  ♦ O’)  (4.38) 

a - b 

as  the  asymptotic  speed  of  the  disturbance  in  a cylindrical  rod.  This 

23 

speed  is  exactly  that  predicted  by  Rayleigh’s  approximate  theory  for 


^Lord  Rayleigh,  Theory  of  Sound  (Cambridge,  1877),  I.  pp.  242-251. 


thin  rods. 

The  simple  expression,  valid  for  all  real  values  of  Cy,  obtained 
for  the  asymptotic  properties  as  {,'L/2Rj  — > , leads  one  to  reexamine 

the  processes  leading  to  equations  (4.26)  and  (4.27)  with  the  idea  of 
modifying  theee  processes  so  as  to  obtain  a segregation  of  the  disturbance 

into  parts  v^hich  will  correspond  exactly  with  the  geometrical  paths  of 
Section  3 regardless  of  the  value  of  ~l)  . 


4.19 

Haring  found  the  equation  (4*27)  which  is  consistent  with  such  a 

segregation  whenever  slnVtr  - 0,  it  is  natural  to  consider  the  possibility 

of  extending  these  results  to  other  values  of  l)  by  interpolation  between 

these  values  cr  some  other  such  process  of  generalisation. 

The  integration  over  s in  equation  (4«27)  is  readily  carried  out  by 

evaluating  residues  with  the  aid  of  Jordan's  lenaa  at  |2V  ♦ l|  poles  of 

order  n^  ♦ n^  * 1.  The  results  so  obtained  will  give  quite  complicated 

combinations  of  exponential  and  algebraic  functions  which  must  then  be 

integrated  over  y.  In  order  to  point  out  the  fruitless  nature  of  any 

effort  at  interpolation,  it  is  only  necessary  tc  consider  the  values  of 

A.  (R  »Z,t~t  ) obtained  from  equation  (4-27)  for  the  initial  rise 
"ln2  ° , 

Just  after  t • t t°  . Since  the  initial  rise  is  controlled  almost 

nln2 

entirely  by  the  Fourier  components  of  high  frequency,  it  is  convenient  to 

expand  T"  (sh,  sk,  R ) in  inverse  powers  of  a and  use  the  leading 
nln2  ° 

term  .xnly  in  getting  an  approximation  to  the  behavior  in  the  neighborhood 
of  the  initial  rise  time.  This  is  easily  accomplished  with  the  asymptotic 
formulae  (4-20)  and  (4*21)  which  are  identical  when  sinV  n « 0. 

The  res\ilt  may  be  written  as 

A <V)  I A t 

A (P.  ,Z,t-t  ) - (2V  ♦ 1)  (-1)  A (R  ,z,t-t  ) (4.39) 

•^n^n^  o ' ' nln2  0 
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t° 

nln2 


< t < t 


♦ t° 

nln2 
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where  £ is  a small  positive  number, 
equation  (3*63). 


and 


A 


-v  (R  ,z,t- 

V>2  ° 


•t  ) 


is  given  by 
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From  a physical  point  of  view,  the  main  feature  of  equation  (4*39) 

l>(nl*n2) 

is  the  phase  changing  factor  (-1)  . There  appears  to  be  no  way 

of  generalization  consistent  with  experiment  in  the  case  y)  • l/2,  which 

has  been  investigated.^  Experimentally  it  appears  that  the  initial  rises 

of  * (R  ,Z,t)  are  of  the  same  sign  as  that  of  the  direct  dilatational 

' 'o  n?'  o 

wave  and  anpear  for  all  values  of  n~  just  as  would  be  expected  of  a plate. 

'iMni«-n2) 

If  one  attempts  to  generalize  by  replacing  (-1)  x by 

cos  £ (n,  + n2)>»n^,  rises  could  only  be  found  for  even  values  of  n^ 

when  n,  - 0 and  7)  - 1/2,  and  these  must  oscillate  in  sign  accordingly 

as  n0  increases  through  such  even  values., 

t)  (ni+n2) 

Since  the  same  factor  (-1)  is  also  contained  in  the  neglected 

tcrm3,  it  is  apparent  that  some  other  approach  is  required  when  sin^n  / 0. 
However,  the  nature  of  thi3  approach  is  not  immediately  obvious. 

The  complementary  solution  for  ■ 1/2  can  be  calculated  from  the 
equations  developed  in  this  section,  but  the  effort  and  time  required 
would  be  prohibitive. 

The  failure  t.o  obtain  a solution  like  that  for  the  plate,  in  which 
individual  terras  of  the  series  correspond  to  each  of  the  various  types 
of  geometrical  paths,  is  certainly  not  to  be  construed  as  clouding  the 
significance  of  such  paths  or  the  existence  of  the  corresponding  wave 
trains,  for  these  wave  trains,  interpreted  according  to  the  principles  of 
geometric  optics,  have  been  used  in  geophysical  prospecting  to  map  all 
kinds  of  curved  surfaces  in  a consistent,  manner  and  in  agreement  with  data 
obtained  from  well  cores.  The  particular  mathematical  met  hod  employed  has 
simply  been  inadequate  to  ace  jap  li  e the  split  tin.-  of  the  overall  !ic 
turbaries  into  senara*  0 o;  ecc'  el’  t he  desired  type.. 
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The  mathematical  methods  ware  adequate  for  the  cases  in  which 
3in)/rr  - 0,  as  equations  of  the  form 


- (1/^L  (-z)€i(z  ' n(  V*  1)/2) 

L.^  (4.40) 

♦ Pt)  (z)^'i(Z  ' U{V  + 1)/2}] 


exist  where 


i l)  - V2 1 

fy<i)  ■ ZH 

m - 0 


1/2  (It,/*  1/2 i - 1/2  ♦ m)l 

(17/  + 1/2 j - 1/2  - m)!  ml  (2iz)m 


(4*41) 


when  is  a positive  or  negative  integer  or  zero.  The  use  of  these  ex- 
pressions in  the  preceding  analysis  gives  equation  (4.27). 

When  7/  is  not  a positive  or  negative  integer  or  zero,  there  appears 
to  be  no  function  p-jy(z)  satisfying  equation  (4. 40)  which  ha3  the  property 
that  the  ratio  p^,  z)  i3  single  valued  and  regular  at  infinity 

(i.e.  a rational  function).  However,  it  is  in  general  possible  to  re- 
present the  Bessel  functions  of  the  first  kind  as  the  limit  of  a sequence 
of  term3  of  the  nature  of  equation  (4-40)  namely 

• . (u,N) 

J-^y.U)  ’ Lim 

N 00 


( 4 . 42 ) 
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whera 


U>  ■ (i/3n;)'-/2  (z/^^rp^.-Oe1'1  ' "("  * 1)/2) 


(4.43) 


♦ P^U)€"i(z  " n(M  + 1)/2) 


u is  an  integer,  and 

N 

nS  (*)  “ nv  * i - n)  7 ] 

>u  n=i 0 

ill  * 2n  «■  1/2)  (u  * n ♦ 1/2)  ( ~l)n 

• p z) 

ni  T(V  + 1 - u - n)  V{V  * n ♦ 3/2)  u 


( 4*  44) 


In  which  the  ratio  (-z)/p“  (z)  i3  obviously  a rational  function. 

This  infinite  process  is  easily  obtained  from  Sonine' s expansion 
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Watson, 


op. 


cit . , 


pp.  139-140. 
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'(0=1 


(m  ♦ 2n  + 1/2)  T\l‘  * n ♦ l/2) 
n\  V('^  * 1 - u - n)  P(  * v.  ♦ 3/2) 


+ 2n+  ^ 


(4.43) 
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which  is  valid  when  u * 1/2,  + l/2,  and  •})  - /i  are  not  negative  inte- 
gers. One  merely  defines  as  the  sum  of  the  first  N terms  of 

the  above  series  with  the  stipulation  that  m is  an  integer  and  applies 
equation  (4. 40)  to  each  term  sc  considered. 

ALthough  these  processes  converge  quite  rapidly,  one  finds  that 

Lim  p*  (2)  does  not  in  general  exist.  Exceptions  occur  when  1) , as 
N — >o°  ^ 

well  a3  u,  is  ?r:  ..ntoger  and  1/  'S-  p..  In  these  cases  the  sequence  termi- 
nates, and  one  find3 


V-> 


P^j  (z) 


1)  zt  M » N > ^ - u 


(,4.40; 


—fl?) 

The  possibility  of  representing  as  the  limit  of  a sequence 

Jl)  M U) 

based  on  the  above  infinite  sequence  suggests  itself.  Defining  S ' ' 

by  making  the  substitutions 


zzc 


_( >\u,fO 


zzc 


I v J(d,N) 


(4.47) 
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i)-y^z) 


I (d-1  ,N ) 


into  equation  (4.1)  and  expanding  the  right-hand  member  in  ascending 
Ro  and  (£  ^RRo,  one  readily  obtains 


powers  of  £ 


ao  00 
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where 


O oo  + i T 


(Ro,Z,t-t  ) 


-i(2V>  1)(1  - 2b2/a2)2(2n)2 

oo  ■+  L ~C 
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nln2 
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, b*h3R 
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^ a3/ 2 a / i w>.  _ _ -is (t-t  -t  ) 

— ; — r — (sh»sk>R  i C 12  dyds  (4.49) 

v 4r3~  1 n.n-  O’ 


> 


van- 


in  v.hich  1”  ^ (sh.sk,  ii  ) is  a rational  function  of  each  of  the 

nln2  o' 

ables  sh,sk.  and  R . It  has  poles  of  order  n,  ■»  n * 1,  considered 

C 1 

as  a function  of  s,  and  the  number  of  these  poles  is  dependent  upon  ~\) , u» 

and  N.  It  is  easily  demonstrated  that  A<y„  ,m,N^(R  ,Z,t-t’ ) - 0,  for 
. nl:  2 0 

I «*»Q  % 

t-t  ■<-  tn  n and  for  any  value  of  'v  , and  that  this  result  is  compatible 
with  equation  (4.27),  when  3inV  tt  * 0,  in  the  sense  that 


^ n.n.  > An,n0  "V^"'  ' 

N — ^ co  1 ■c  12 


(4.50) 


as  a result  of  equation  (4*46). 

If  this  process  converges  when  sin  n f 0,  equation  (4-50)  would  be 

a much  more  satisfactory  definition  of  (R  ,Z,t-t  ) than  that  pro- 

/'n^n^  o' 

vided  by  equation  (4  2b) , as  this  impulse  response  would,  then  always  corre- 
spond to  the  appropriate  geometrical  path. 
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V 

A 

Although  the  integration  over  3 in  equation  (4-49)  ccuid  be  carrxed 
out  formally  in  terms  of  residues,  the  number  and  location  of  the  poles 
depends  upon  N in  such  a way  that  the  author  has  been  unable  to  evaluate 
the  limit  involved  in  equation  (4.50)  or  prove  that  it  exists.  The  physical 
role  of  the  geometric  paths  leads  one  to  expect  the  limit  to  exist,  but 
there  is  very  little  hope  of  reducing  the  calculation  to  processes  of  suf- 
ficient simplicity  to  warrant  further  consideration  in  this  paper. 


5 . CONCLUSION 


A theoretical  investigation  has  been  made  of  the  propagation  of 
elastic  pulses  through  roda  and.  plates.  This  investigation  is  by  no 
means  exhaustive  as  it  has  been  limited  to  a plane  longitudinal  drive. 
Transverse  drives  and  point-source  drives  are  equally  interesting  from 
the  point  of  view  of  making  physical  measurements  of  elastic  constants, 
at*  are  the  parallel  problems  involving  nor.-isotropic  media. 

In  the  case  of  the  plate,  these  mathematical  methods  are  capable  of 
treating  the  transverse  drive,  point-source  drives,  and  the  non-isotropic 
media  with  minor  alterations.  The  author  intends  to  treat  these  and  certain 
related  geophysical  problems  in  future  papers. 

In  the  c*a3e  of  the  rod,  it  is  clear  that  some  other  mathematical  tool 
must  be  annl ied  to  the  much  more  complicated  functions  wi  ich  »re  involved. 
The  primary  objective  of  this  research,  a quantitative  wave  treatment  of 
the  experiments  reported  by  Hughes,  Pondrom,  and  Mims,  has  not  been  com- 
pletely attained.  At  present  the  best  approach  is  to  consider  the  effects 
that  appear  in  connection  with  rods  in  terms  of  an  analogy  with  those  that 
appear  in  connection  with  the  plate.  This  is,  at  best,  only  a qualitative 
analysis. 

For  a longitudinal  drive  on  a plate,  it  has  been  shown  that  the 
boundaries  of  the  plate  produce  reflected  or  echo  wave  trains  which  corre- 
spond to  each  of  the  paths  predicted  by  the  methods  of  geometric  optics. 

It  is  found  that  the  critical  anglo  paths  should  be  the  most  evident 
experimentally  as  the  associated  wave  trains  have  a much  more  abrupt  start. 
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It  ia  further  found  that  these  wave  trains  interfere  with  one  another  in 
such  a way  that,  as  the  ratio  of  the  distance  of  transmission  to  the 
thickness  of  the  plate  becnjies  infinite,  the  entire  disturbance  becomes 
asymptotically  a disturbance  traveling  with  the  velocity  cq  which  is 
classically  the  longitudinal  velocity  of  propagation  for  a plate  of  zero 
thickness.  These  results  are  obtained  independent  of  any  assumptions  re- 
garding the  frequency  spectrum  of  the  pulse. 

Oy  analogy,  one  should  expect  longitudinal  drive  on  a rod  to  result 
in  wave  trains  which  correspond  to  each  of  the  paths  preiicted  by  means 
of  geometrical  optics,  but  that  those  wave  trains  corresponding  to  the 
critical  angle  paths  and  the  direct  dilatational  path  will  be  the  mo3t 
easily  detected.  As  the  ratio  of  the  distance  of  transEiission  to  the 
diameter  becomes  infinite,  interference  of  the  above  wave  trains  should 
produce  the  asymptotic  result  of  a single  wave  train  propagated  with  the 
speed,  found  by  Rayleigh  for  a rod  of  zero  radius.  It  is  also  to 

be  expected  that  the  frequency  spectrum,  or  shape  of  the  driving  pulse, 
is  of  no  importance  to  the  theory.  All  of  the  conclusive  findings  of 
Section  are  in  accord  with  this  analogy. 

The  experimental  observation  of  the  various  wave  trains  is  complicated 
by  the  multitude  of  ways  in  which  interference  can  take  place.  This  be- 
comes particularly  difficult  when  the  duration  of  the  driving  pulse  F(t) 
is  greater  than  the  differences  between  the  minimum  group-transit  times 
of  several  of  the  various  reflected  wave  trains.  This  condition  existed 
in  the  experimental  work  reported  by  Hughes,  Pondrom,  and  Mlms^  and  ef- 
fectively prevented  tr.o  identification  of  any  wave  traina  except  those 
corresponding  to  the  direct  dilatational  and  critical  angle  paths  which 
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correspond  to  the  ccntoura  n . The  more  general  paths  corresponding 

to  f with  n,  small  but  not  zero  generally  give  rise  to  slow  starting 

transients  as  pointed  out  in  Section  3*  However,  as  (Z/2Ro)  is  made  larger, 

these  rises  become  more  near?.y  like  those  corresponding  to  the  critical 

angle  paths  and  should  then  he  mor  ' readily  detected.  Unfortunately,  the 

minimum  group  transit  times  t°  — ► t°  as  (Z/2R  ) — ► °°  , and  the 

& ' nln2  0 n2  o' 

wave  trains  corresponding  to  the  more  general  paths  interfere  with  other 

wave  trains  under  the  same  conditions  that  they  give  fairly  sharp  rises 

unless  F(t)  is  sufficiently  short  in  duration  to  avoid  this  interference. 

This  point  is  illustrated  by  Fig.  5*1  which  is  a reproduction  of  the 

oscillographic  traces  for  a cold-rolled  steel  rod  3 * 61,5  cm.  in  diameter 

and  5 -OP  cm.  in  length.  For  this  sample  of  steel,  a was  found  to  be  5880 

raetera/sec. , b was  found  to  be  3203*5  meters/sec.,  and  the  value  .289  was 

obtained  for  CT , 

The  primary  dilatational  disturbance  begins  at  point  A,  8.69  Msec, 
after  the  start  of  the  driving  pulse,  which  had  an  appreciable  amplitude 
for  some  U Msec.  Since  the  minimum  group-transit  rimes  t^  q associated 

with  the  contours  are  given  by 

“1° 


t° 
nn  o 


(1/a) 
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x O 
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as  the  entire  path  is  traversed  as  a dilatntiorml  wave,  we  thua  find 


>45  cm.  in  diameter  and  5 *08  cm > in  length 
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Confining  our  attention  to  t, 


we  find 


1 o 


/ 8.69 


(3. 645/5- 08) 
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or 


t:  « 10.696  Msec. 

± o 


Thuej  the  corresponding  train  of  wayes  would  interfere  with  the  direct 
dilati tional  wave  train  over  half  of  its  length.  Ao  n^  increases,  scune 
of  the  corresponding  wave  trains  would  be  expected  to  be  resolved  in  terms 
of  time  but  would  not  be  expected  to  have  a sharp  enough  rise  to  be  de- 
tect. cu. 

The  obvious  arrival  at  point  3 is  due  to  a critical  angle  path  corre- 
sponding to  the  contour  1^  1 . The  arrival  time  t^  ^ is  found  to  be  13.61 
Msec,  after  tne  s ..rt  of  the  driving  pul3e  arid  is  in  agreement  with  calcu- 
lations based  upon  equation  (3*52).  The  wave  train  corresponding  to  the 
minimum  group-transit  time  t-,  , which  is  21.23  uaoc  as  calculated  from 

-k  J. 


r r 
?•-> 

equations  (3*46) , (3>U7)>  and  (3*40),  and  thus  it  should  interfere  with 
the  wave  train  corresponding  to  t^ 

These  conditions  become  much  aore  aggravated  as  (z/2Kq)  increases. 
Since  th ; remainder  of  the  rods  considered  have  larger  length  to  diameter 
ratios  than  the  one  considered  Above,  it,  is  apparent  that  there  was  little 
likelihood  of  identifying  any  wave  trains  except  those  corresponding  to 
the  direct  dilatational  and  critical  angle  paths  which  always  have  sharply 
rising  initials. 

A more  promising  set  of  circ;,p  stances  for  the  observation  cf  the 

more  general  wave  trains  corresponding  to  the  minimum  group-transit  times 

t°  with  n,  4 0 can  be  obtained  by  either  shortening  the  duration  of 
nln2  1 

the  driving  pulse  or  by  increasing  both  the  length  and  the  diameter  of 

the  rods  to  be  considered  by  a factor  of  approximately  ten  so  as  to  get  a 

—o 

greater  separation  of  neighboring  values  of  t 

nln2 

The  ’’attar  alternative  is  not.  very  acceptable  to  this  laboratory  as 
the  fundamental  reason  for  developing  this  method  of  measuring  velocities 
of  elastic  waves  was  to  get  a system  in  which  the  rod-shaped  sample  and 
the  necessary  driving  and  detecting  crystals  could  be  placed  in  a small 
volume  inside  of  a high-pres3ure  chamber.  Thus,  with  one  aet-up,  one- 
can  measure  the  velocities  cf  rotational  and  dilatational  waves  for  various 
hydrostatic  pressures  and  temperatures.  It  is  consequently  desirable  to 
pursue  the  former  alternative,  and  efforts  are  now  being  made  to  reduce 
the  duration  of  F(t)  and  increase  its  amplitude  so  as  to  increase  the 
number  of  arrivals  capable  of  being  identified. 


,3^ 
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The  question  of  precisely  what  is  being  measured  in  these  experiment# 

is  closely  tied  in  with  the  nature  of  the  driving  pulse.  If  F(t)  - 0 for 

t < 0 and  is  non-zero  for  any  short  time  interval,  it  is  a well  recognized 

fact  that  this  function  must  be  considered  to  have  a continuous  distribution 

of  frequencies  running  from  0 to  «o , All  elastic  materials,  for  one  or  more 

reasons,  are  expected  to  become  dispersive  as  the  frequency  of  a simple 

hartaonic  disturbance  becomes  sufficiently  large.  Thus,  the  velocities  a 

and  b are  complex  functions  of  the  frequency  which  is  represented  in  the 

analysis  of  this  paper  by  the  symbol  e.  The  velocities  measured  area 

apparently  those  associated  with  the  first  arrival  of  an  abruptly  initiated 

disturbance  or  ,:wave-front''  velocities.  Paralleling  the  reasoning  of 

9 

Sommerfeld  and  Brillouin  , one  might  say  that  the  velocities  measured  are 
really  the  limits  of  a and  b as  s — > oo  . Such  a statement  is  rather  naive 
since  the  results  of  any  such  measurement  are  very  likely  to  be  a function 
of  the  sensitivity  of  the  detecting  mechanism.  Without  a complete  and  ex- 
perimentally acceptable  theory  of  dispersion,  one  cannot  say  when  the  de- 
tector is  sensitive  enough  to  detect  the  very  first  arrival.  There  is  some 
reason  to  expect  that  an  extremely  small  part  of  the  energy  of  an  abruptly 
starting  elastic  disturbance  is  propagated  with  the  3peed  of  light*, 

In  the  absence  of  a complete  theory  of  elastic  dispersion,  no  exact 
description  of  the  measured  velocities  can  be  given.  One  can  only  state 
that  the  measured  velocities  are  the  apparent  "wave-front/'  velocities  as 
observed  with  apparatus  of  specified  characteristics.  It  is  to  be  expect- ed 
that  an  acceptable  theory  of  dispersion  will  result  from  the  study  of  the 
re;..i3t.s  of  ench  experiments  as  these. 
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